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NESTED  270  DEGREE  ELECTROSTATIC  ANALYZER 


Initial  design  studies  and  computer  modeling  calcu¬ 
lations  were  started  on  an  electrostatic  analyzer  (ESA) 
system  using  nested  270  degree  quadr i spheres .  The  outer  set 
of  spherical  deflection  plates  will  be  used  for  ion  analysis 
and  the  inner  pair  for  electron  measurements.  Initial  design 
assumptions  have  been  that  the  ion  sensor  will  require  a 
sensitivity  of  lE-03  cm  squared  steradian  and  the  electron 
sensor  will  be  approximately  ten  to  one  hundred  times  less 
sensitive  than  the  ion  sensor.  This  sort  of  system  could  be 
used  in  a  low  to  medium  earth  orbit  where  the  ion  densities 
are  ten  to  a  hundred  times  less  than  the  ambient  electron 
flux.  This  sensitivity  differential  would  yield  similar 
ion/electron  counting  statistics  for  that  environment. 

The  dimensions  and  characteristics  of  the  deflection 
system  that  was  modeled  are: 


Ion  Plates 


Electron  Plates 


Outer  Radius 
Inner  Radius 
Aperture  = 

Outer  Radius 
Inner  Radius 
Aperture  = 


=  2.3  inches 
=  1.7 
0.067  " 

=  1.125  inches 

=0.875 

0.039 


A  preliminary  estimate  of  the  response  of  this  system 
indicates  that  the  ion  plates  could  sample  an  incoming  flux 
within  an  angular  acceptance  of  120  degrees  by  12  degrees. 

If  the  broad  angle  response  were  divided  into  ten  12  degree 
sectors,  then  these  individual  sectors  would  have  a 
sensitivity  of  lE-03  cm  squared  steradian,  the  initial  design 
goal.  The  plates  would  have  an  analyzer  constant  of  3.33  and 
an  energy  resolution  of  1.7%. 


The  electron  plates  would  likewise  view  more  than  120 
degrees,  but  the  narrow  fan  angle  would  be  10.5  degrees. 
Assuming  an  identical  12  degree  sector  for  the  electrons 
yields  a  sensitivity  of  3E-04,  an  analyzer  constant  of  4  and 
an  energy  resolution  of  2%.  The  geometric  factor  of  the 
electron  deflection  plate  pair  can  be  adjusted  by  varying  the 
aperture  diameter.  An  inlet  aperture  of  .0039  inches  reduces 
the  response  sensitivity  to  2.6E-06. 

These  initial  estimates  were  manually  calculated  and  do 
not  take  into  consideration  detector  efficiencies,  angular 
modulation  of  the  effective  aperture,  deflection  plate  field 
inhomogeneities,  and  other  second  order  effects  on  real  ESA 
systems.  Tnese  preliminary  plate  and  aperture  dimensions  are 
being  run  on  a  computer  model  that  contains  most  of  these 
variables  and  various  performance  and  design  trade  offs  are 
being  considered. 
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The  nested  spherical  design  offers  a  compact  measurement 
module,  but  the  mounting  and  electrical  isolation  of  the 
deflection  plates  and  electron  multipliers  presents  a 
mechanical  design  challenge.  Several  plate  mounting  methods 
have  been  attempted  and  a  method  that  provides  both  the 
isolation  and  mechanical  rigidity  necessary  for  space  flight 
has  been  identified. 

In  order  to  extend  the  angular  coverage  of  this  sensor, 
a  system  using  motor  driven  rotary  platforms  has  been 
investigated.  Two  rotating  nested  pairs  will  achieve  full  2 
pi  response. 

This  sensor  design  will  produce  large  quantities  of  data 
that  could  require  on  board  data  processing  if  there  is  not  a 
substantial  telemetry  bandwidth  to  transfer  the  raw  data  to  a 
ground  station.  A  sensor  control  device  is  being  designed  to 
increase  the  versatility  of  the  ESA.  The  work  done  on  this 
control  unit  is  covered  elsewhere  in  this  report. 

Since  this  ESA  configuration  can  measure  simultaneous 
ion/electron  fluxes  made  over  a  broad  magnetic  pitch  angle 
distribution,  the  experiment  offers  a  unique  opportunity  to 
examine  the  raw  data  stream  for  statistically  significant 
particle  bunching  that  may  occur  in  certain  wave/particle 
interactions.  To  make  this  type  of  measurement  the 
feasibility  of  including  a  particle  correlator  in  the  data 
processing  has  been  studied.  A  University  of  Sussex  expert 
in  this  field  (Dr.  Paul  Gough)  has  been  contacted  about  this 
possibility. 

Vendors  of  several  types  of  electron  multipliers  have 
been  contacted  to  identify  suitable  candidates  for  the  sensor 
system  multiplier.  Laboratory  experiments  to  evaluate  some  of 
these  microchannel  electron  multiplier  arrays  have  been 
conducted.  We  are  looking  for  a  source  of  large  area 
multiplier  plates  that  have  good  output  pulse  height 
distribution  and  are  rugged  enough  for  a  space  experiment. 
Samples  from  two  manufacturers  have  been  obtained  and  are 
being  evaluated. 

The  possibility  of  including  a  data  recording  device  as 
part  of  the  experiment  is  also  being  explored.  Tape 
recorders,  Winchester  type  hard  disks,  video  recorders,  and 
optical  disk  recorders  are  being  researched  to  identify 
manufacturers  of  ruggedized  systems  capable  of  space  flight. 
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ION  APERTURE  ELECTRON  APERTURE 


4 


NESTED  PAIR  ESA 


DATA  PROCESSING  UNIT  FOR  270  DEGREE  NESTED  SENSORS 


The  Data  Processing  Unit  (DPU)  is  in  the  conceptual 
design  stage.  Because  the  data  requirements  remain  unknown, 
the  present  configuration  is  a  dual-CPU  system.  This 
con  figuration  should  al low  sufficient  processing  overhead  for 
the  estimated  data  requirements. 

In  this  configuration,  one  CPU  is  dedicated  to  data 
acquisition  from  the  sensors  and  particle  correlator,  and 
controls  the  high  voltage.  It  performs  all  of  the  data 
analysis  and  transfers  only  the  relevant  data  to  the  second 
CPU. 

The  second  CPU  is  responsible  for  formatting  the  data 
received  from  the  first  CPU  and  transferring  it  to  the  data 
recorder.  It  is  also  responsible  for  controlling  the  motion 
system  and  monitoring  the  environment.  Communications  with 
the  spacecraft  are  controlled  through  this  CPU. 

Each  CPU  has  its  own  local  program  and  data  memory,  and 
other  necessary  peripherals.  The  only  interconnection  between 
the  two  CPUs  (other  than  power)  is  a  16-bit  bidirectional 
data  bus.  This  bus  allows  the  transfer  of  timing  information, 
mode  control,  and  sensor  data  between  the  two  CPUs. 
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DATA  PROCESSING  UNIT  FOR  NESTED  SENSORS 


FIGURE  QRl- 


FIGURE  QRl- 


High-Voltage  Optocoupler 


A  high-voltage  optocoupler  is  being  developed  for  use 
within  an  instrument's  high-voltage  power  supply.  The 
optocoupler  will  perform  as  a  current  controlled  current 
switch.  Typical  applications  include  regulating  the  potential 
between  deflection  plates  or  across  a  microchannel  plate. 

The  output  switch  can  be  modeled  as  a  5000  volt 
photodiode.  The  input  is  a  string  of  light  emitting  diodes 
(LEDs) .  Infra-red  light  generated  by  the  LEDs  causes  the 
reverse  biased  photodiode  to  'leak'  current  by  the 
photoelectric  effect.  The  output  current  ('leakage')  can  be 
controlled  by  the  input  current  to  the  LEDs. 

In  practice,  a  5000  volt,  low-leakage  photodiode  does 
not  exist.  A  stack  of  1000  volt  transistors  or  diodes  is 
being  proposed  as  an  equivalent.  Custom  LEDs  are  being 
examined  to  improve  light  generating  efficiency.  The  unit 
will  be  assembled  into  a  form  of  hybrid  circuit. 

Various  optical  couplings  are  also  being  developed.  The 
couplings  must  be  able  to  gather  and  transmit  a  maximum  amount 
of  light  as  well  as  withstand  at  least  a  10,000  volt  isolation 
voltage.  Parabolic  reflectors,  light  pipes,  fiber  optic 
faceplates,  and  several  silicone  encapsulations  are  presently 
under  test  as  such  couplings. 

Some  goals  of  this  development  are  to  obtain  an 
optocoupler  that  can  withstand  an  input-output  isolation 
voltage  of  10,000  volts,  hold  an  output  differential  voltage 
of  5000  volts,  and  pass  500  microamperes  of  output  current 
with  only  20  milliamperes  of  input  drive  current. 

The  following  two  figures  show  the  schematic  and  layout  of 
the  latest  design  approach.  This  layout  requires  that  all 
parts  are  mounted  on  a  single  piece  of  alumina.  A  silicone 
'light  pipe'  is  used  for  coupling.  The  optical  path  is  formed 
by  the  curved  shape  the  silicone  forms  due  to  surface  tension. 
An  outer  coating  of  white  alumina  filled  silicone  also  aids 
internal  reflection. 

This  latest  configuration  has  a  light  collection 
efficiency  of  200  microamperes  output  for  20  milliamperes 
input.  Output  blocking  voltage  is  initially  5000  volts,  but 
decays  after  several  hours.  Understanding  the  reasons  for 
this  decay  has  become  a  top  priority. 

The  demanding  breakdown  and  isolation  voltages  will  make  this 
part  unique,  unlike  anything  available  commercially.  High 
speed  and  compact  size  are  two  predominate  features  that  will 
improve  existing  high-voltage  power  supply  designs. 
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FIGURE  QRl 
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SENSOR  DESIGN 


The  development  of  a  motorized  nested  hemispherical 
electrostatic  analyzer  results  in  a  sensor  system  that  can 
view  more  than  two  pi  solid  angle.  Because  of  this  enormous 
viewing  angle,  it  is  difficult  to  mount  the  sensor  on  a 
spacecraft  such  that  parts  of  its  acceptance  angle  will  not 
be  obscured  with  spacecraft  structural  members  or  adjacent 
experiments.  Surfaces  near  or  in  the  sensors  field  of  view 
can  effect  the  data  either  by  shadowing  or  by  charge 
accumulations  that  generate  electric  fields  that  will 
deflect  incoming  particles. 

Amptek  has  been  using  computer  simulations  and  graphics 
methods  of  modeling  the  sensor  system  and  its  associated 
fields  of  view.  These  methods  allow  the  three  dimensional 
modeling  of  surfaces  and  structures  in  the  vicinity  of  the 
sensor  and  the  identification  of  the  degree  of  obscuration. 
This  method  simplifies  optimizing  the  sensor  viewing,  the 
range  and  rate  of  mechanical  sweeping,  and  the  necessary  data 
processing  algorithms  needed  to  interpret  data  sets  taken 
while  the  view  is  obscured. 

Satellite  structures  and  sensor  mounting  on  them  can 
often  accommodate  large  viewing  angles,  but  experiments  flown 
in  the  STS  shuttle  bay  will  have  limited  ranges  of  mounting 
possibilities.  Because  the  shuttle  bay  seems  to  be  a  worst 
case  example  of  limited  and  obscured  viewing,  the  modeling 
technique  was  used  to  assess  a  pair  of  nested  hemisphere 
ESA's  mounted  in  the  shuttle  bay.  Figures  1  and  2  show 
examples  of  some  of  these  viewing  problems.  An  experiment 
mounted  on  a  pallet  that  is  above  the  bay  door  sill  will  have 
reasonably  good  viewing.  The  principle  obscuration  would 
result  from  adjacent  experiments  and  pallets. 

Amptek  has  been  discussing  microchannel  plate  designs 
and  configurations  with  a  manufacturer.  Figure  3  shows  a 
possible  arrangement  which  will  allow  a  single  shape  to  map 
both  the  electron  and  ion  sensor  responses.  If  a  single  shape 
can  be  developed,  the  cost  of  manufacturing  the  plates  should 
be  reduced  and  the  necessary  inventory  of  spare  plates  needed 
for  future  flight  projects  would  be  minimized.  It  is 
important  to  make  an  early  definition  the  microchannel  plates 
because  they  are  state-of-the-art  developments  that  have  long 
lead  times  associated  with  their  manufacture. 
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Sketch  of  how  this 

shape  can  cover  both  ions  and  electrons 


Large  area  micro-  ___ 
channel  plate  shape. 
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HIGH  VOLTAGE  OPTOCOUPLER 


The  Optocoupler  project  is  aimed  at  developing  a  small  and 
efficient  method  of  controlling  the  output  stage  of  a  high- 
voltage  (5-lOkV)  source.  Typical  applications  include  biasing 
deflection  plates. 

The  search  for  an  LED  emitter  for  the  input  stage  has  ended 
with  a  custom  design.  Applications  for  the  Optocoupler  will 
require  significant  light  output  from  the  LEDs  for  several 
years.  Typical  devices  are  not  very  efficient  and  would  have 
deteriorated  rapidly  when  driven  to  provide  the  necessary 
output  power.  The  custom  approach  started  by  choosing  a 
manufacturer  with  an  existing  high  efficiency  part. 
Modifications  were  made  to  a  standard  LED  chip  to  quadruple 
its  existing  surface  area.  This  effectively  cut  the  strain  in 
each  junction  by  75%.  Degradation  appears  to  react 
exponentially  with  strain,  so  the  new  LEDs  are  expected  to 
last  well  over  10  years  under  moderate  service.  Prototypes 
have  been  delivered  from  the  manufacturer.  The  custom  LED 
chip  is  depicted  in  figure  HV-1. 


Numerous  attempts  have  been  made  to  obtain  a  stable,  high 
voltage,  high  gain  silicon  junction  for  the  output  stage. 
Studies  have  been  conducted  with  both  diodes  and  transistors. 
The  high  sensitivity  to  contamination  of  high  voltage  devices 
has  posed  many  problems.  In  the  interest  of  obtaining  an 
extremely  robust  optocoupler,  a  custom  diode  is  also  being 
developed.  This  diode  design  is  heavily  passivated  with 
glass.  It's  mesa  structure  aids  in  the  passivation  technique 
and  helps  prevent  premature  'hotspots'.  The  large  open  area 
on  the  top  of  the  chip  is  optimized  for  IR  light  collection. 
Breakdown  voltages  for  these  devices  should  range  from  600V 
to  800V,  with  a  dark  current  leakage  of  less  than  20nA. 
Delivery  of  the  first  prototypes  is  expected  at  the  end  of 
May,  1988.  The  high  voltage  diode  chip  is  shown  in  figure 
HV-2 . 


The  process  of  optically  coupling  the  IR  LEDs  to  the  light 
sensitive  region  of  the  high  voltage  diodes  can  be  approached 
in  many  ways.  The  most  promising  technique  appears  to  be  with 
the  introduction  of  a  fiberoptic  faceplate.  Figure  HV-3 
roughly  describes  a  method  for  mounting  all  the  critical 
parts  together.  Both  the  LEDs  and  the  diodes  are  mounted  on 
separate  substrates  and  isolated  by  the  fiberoptic  glass. 
Fiberoptics  are  used  to  transmit  the  'image'  of  the  LEDs  onto 
the  diodes  with  a  minimum  amount  of  light  spreading  or  loss. 
This  design  couples  a  maximum  amount  of  light  between  the 
active  devices  while  providing  high  voltage  isolation.  The 


entire  assembly  is  potted  in  a 
for  hermeticity  and  dielectric 
been  omitted  for  clarity. 

These  processes  will  produce  a 
device,  capable  of  efficiently 
frequencies . 


high  purity  silicone  elastomer 
strength.  Connections  have 

very  compact  regulating 
controling  currents  at  high 
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Ooto 


Amptek  par't  numbers : 


CG3. 0/6-30:  Dark  current  <  30nA  ?  600V 
CG3. 0/8-30:  Dark  current  <  30nA  §  300V 
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CASSETTE  TAPE  SUBSYSTEM 


A  high  capacity  tape  drive  is  being  developed  to  address  the 
problem  of  acquiring  massive  volumes  of  data.  Modern 
recoverable  experiments,  such  as  those  aboard  shuttle,  have 
the  opportunity  and  technology  to  generate  enormous  amounts 
of  data.  Particularly,  imaging  data.  Down  links  via  radio  are 
often  insufficient  for  the  amount  of  information  available. 
The  Exabyte  EXB8200  Cassette  Tape  Subsystem  (CTS)  may  provide 
an  economical  solution  for  storing  up  to  2  giga  bytes 
(2,000,000,000)  of  data. 

The  basic  tape  transport  mechanism  uses  helical  scanning  to 
achieve  its  high  density.  This  technique  has  been  used  for 
many  years  in  analog  designs  such  as  video  tape  recorders 
(VCR).  The  major  breakthrough  with  the  CTS  is  the  ability  to 
store  digital  data.  In  theory,  every  dot  of  every  picture 
produced  by  a  VCR  could  be  a  digital  bit.  Picture  dropouts 
and  other  tape  errors,  however,  are  always  present  and 
require  special  consideration.  Extensive  error  correcting 
processes  must  be  employed  and  are  accomplished  with  several 
application  specific  integrated  circuits  (ASIC).  The  recent 
availability  of  large  scale  ASICs  has  made  this  system 
feasible  for  reliable  digital  data  storage. 

Presently,  the  CTS  is  being  evaluated  for  its  electrical  and 
mechanical  integrity. 

Electrical  tests  will  determine  the  reliability  of  the  stored 
data  and  thus  suggest  appropriate  applications.  Initial 
results  indicate  that  the  CTS  can  at  least  hold  imaging  data, 
where  an  occasional  bit  error  would  go  relatively  unnoticed. 

A  more  stringent  requirement  would  be  to  store  processor 
instruction  codes.  The  feasibility  of  this  application  has 
yet  to  be  determined. 

Mechanically,  the  tape  transport  mechanism  is  quite 
complex.  Tests  of  it's  response  to  a  vibration  input  have 
been  planned.  Barry  Controls  Corp.  has  also  been  contacted  to 
help  with  a  vibration  mount,  if  required.  Enclosing  the 
entire  unit  within  a  hermetic  enclosure  is  planned. 

If  tests  prove  successful  and  a  suitable  enclosure  and  data 
link  can  be  provided,  the  CTS  will  be  a  most  desirable 
storage  warehouse  for  recoverable  experiments. 
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E.'  G  8200  CTS  Product  Specification 


EXABYTE 


Product  Description 

I  he  EXB-82(X)  CTS,  shown  in  Figiirc()R 2-7  consists  of  an  8mm  tape  transport 
mechanism  and  recording  channel,  servo,  formatter,  controller,  interface 
electronics,  software,  and  package  parts  designed  and  produced  by  EXABYTE. 
I  fie  product  is  a  true  digitaJ  data  storage  device,  derived  from  8mm  video  record¬ 
ing  technology,  with  performance  improvements  and  many  additional  functions 
nc ccss.'iry  for  data  processing  purposes. 


FIGURE  QR2-7EXn-8200  Cartridge  Tape  Subsystem 
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Sensor  Design 


The  feasibility  of  a  nested  tri-quadraspheric 
electrostatic  analyzer  continues  to  be  studied.  An  entire 
instrument  complement  necessary  to  control  and  archive 
measurements  is  being  considered.  This  approach  is 
necessary  in  order  to  utilize  the  vast  scope  of  possible 
measurements  that  may  be  made  by  this  type  of  detector 
array.  The  accompanying  drawings  show  the  current  state  of 
these  support  components  and  their  possible  use  on  a  shuttle 
mission . 

Preliminary  thermal  analysis  of  the  instrument 
compliment  indicate  the  need  for  thermal  blankets  and 
possibly  heaters  to  maintain  the  orbital  thermal  excursions 
to  acceptable  levels  for  the  on  board  electronics. 
Preliminary  vibration  analysis  is  being  used  to  further  the 
sensor  and  support  components  design.  Particular  attention 
is  being  paid  to  the  mounting  and  electrical  isolation  of 
the  spherical  deflection  plates  and  the  delicate  micro- 
channel  plate  electron  multipliers. 

Work  has  begun  on  upgrading  the  vacuum  chamber  for  the 
testing  of  the  nested  hemispheres.  A  thermal  shroud  has  been 
designed,  ordered,  and  is  being  installed  to  improve  the 
space  simulation  capabilities  of  the  chamber.  A  small 
electron  gun  has  been  built  and  tested.  It  operated  from  a 
ultra-violet  light  source  used  to  illuminate  an  aluminum 
photo-cathode.  The  resulting  secondary  electrons  were 
accelerated  to  form  a  mono-energetic  electron  beam.  The 
systems  worked  quite  well,  but  failed  when  operated 
continuously  for  a  few  days.  The  UV  lamp  overheated  in  the 
vacuum  and  generated  large  quantities  of  outgassing  and 
contamination.  Alternate  electron  sources  as  well  as  ion 
sources  are  being  considered.  Methods  of  beam  density 
modulation  are  also  being  studied  to  allow  dynamic  testing 
of  a  particle  correlator  design. 
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FIGURE  QR3-5 
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Micro -channel 
electron  multiplier 
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FIGURE  QR3-8 
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DATA  PROCESSING  UNIT  FOR  270  DEGREE  NESTED  SENSORS 


Part  selection  for  the  Data  Processing  Unit  is  nearly 
complete.  The  80C86  microprocessor  family  from  Harris  is 
being  used,  along  with  HM65642  and  HM6617  memory  devices, 
also  from  Harris.  The  D400  from  Amptek  is  being  used  for  the 
IG-bit  accumulators,  and  an  analog-to-digital  convertor 
(ADC)  from  PMI  will  be  used  for  housekeeping  information. 

All  logic  integrated  circuits  will  come  from  Texas 
Instruments . 

The  design  of  the  processor  boards  and  the  accumulator 
board  is  also  nearly  complete.  After  the  operation  of  the 
memory  decoding  system  is  verified,  the  layout  of  printed 
circuit  boards  will  begin. 

Part  selection  is  now  under  consideration  for  the 
Particle  Correlator.  Several  microprocessor  families  are 
under  consideration,  including  the  Inmos  Transputer,  the 
Harris  80C286  with  an  80287,  and  the  combination  of  a  Harris 
80C86  with  a  high-speed  math  coprocessor. 

DPU  components  which  have  not  yet  been  designed  include 
the  particle  correlator,  the  motion  control  interface,  the 
high-voltage  control  interface,  the  spacecraft 
communications  interface,  and  the  power  supply. 
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SPREE  CPU  BLOCK  DIAGRAM 


Project  NO.F19628-87-C-0094  Date:4-28-88  jOwg.  No .  D:  \SPREE\QTR  ID 


SPREE  FLIGHT  RECORDER 


I.  Introduction 

The  proposed  Shuttle  Potential  and  Return  Electron  Experiment 
(SPREE)  will  generate  enough  data  to  mandate  a  flight  recorder  in 
addition  to  a  radio  down  link.  The  many  features  of  SPREE  will 
require  data  storage  in  excess  of  a  gigabyte  (1000  megabytes). 

This  large  amount  of  Information  can  not  be  transmitted  with  the 
limited  bandwidth  generally  provided.  Therefore,  non-critical 
information  will  be  available  only  from  the  flight  recorder. 

Power  constraints  prohibit  the  use  of  several  low  capacity 
recorders.  Instead,  Amptek  is  adapting  new  high  density  recording 
media,  thus  requiring  only  one  or  two  recording  devices.  Media 
include  helical  scan  tape  recorders  (similar  to  VHS  video  rec 
orders)  and  optical  'Write  Once,  Read  Many'  (WORM)  disks. 

II.  Background 

A.  Optical  Disks 

One  of  the  newest  media  technologies  is  the  polyester  optical 
disk,  better  known  as  'laser  disk'.  Current  models  can  write  to 
the  disk  once  by  burning  digital  'spots',  then  read  back  the  data 
nearly  an  infinite  number  of  times.  The  data  is  not  prone  to 
magnetic  erasure  and  has  an  archive  life  exceeding  ten  years. 
Future  designs  will  be  able  to  write  and  erase  data  several  times 
by  means  of  localized  thermal  reflow  of  the  optic  media. 

Mountain  Optotech  Incorporated  manufactures  a  ruggedized  200 
megabyte  optical  disk  drive  for  $6749.  Survivability 
specifications  for  this  unit  include  vibration  tolerance  of  13g 
rms  between  20  and  2000  Hz,  and  shock  to  30g  for  11  msec.  An 
extended  performance  version  intended  for  spacecraft  applications 
expands  the  vibration  specification  to  22g  rms  over  the  same 
frequencies  and  a  has  temperature  range  of  -55  to  90  degrees 
Celsius . 

Despite  its  rugged  specifications,  the  200  megabyte  capacity  does 
not  make  the  Mountain  Optotech  drive  a  first  choice  for  a  flight 
recorder . 

B.  Exabyte  Cartridge  Tape  Subsystem  (CTS) 

Exabyte  has  combined  an  8mm  Sony  camcorder  mechanism  with 
proprietary  Application  Specific  Integrated  Circuits  (ASICs).  The 
result  is  a  2  gigabyte  tape  drive  that  fits  in  a  standard  full 
height  PC  drive  slot  and  consumes  less  than  20  watts  of  power. 

Data  is  transferred  via  a  Small  Computer  Systems  Interface  (SCSI) 
port  at  a  sustained  rate  of  1.8  megabits  per  second. 

Storage  density  is  determined  by  the  techniques  used  in  reading 
and  writing  to  the  tape.  The  size  of  a  data  bit's  magnetic 
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signature  on  the  tape  can  not  be  smaller  than  the  air  gap  in  the 
read  head.  Earlier  data  recorders  using  the  nine  track  technique 
pull  the  tape  past  stationary  heads.  The  maximum  data  rate  can  be 
determined  by  the  number  of  finite  magnetic  transitions  that  can 
occur  in  an  inch  of  tape  multiplied  by  the  tape  speed  in  inches 
per  second  ( ips ) .  With  helical  scanning,  the  head  assembly  spins 
diagonally  over  a  moving  tape,  creating  an  illusion  to  the  head  of 
a  much  faster  tape  speed.  This  process  produces  a  micro  stepping 
of  adjacent  helical  tracks,  yielding  a  very  high  density  and 
allowing  the  tape  to  move  slower. 

Price  of  commercial  units  is  $3,500. 

III.  Development 

Amptek's  goal  is  to  produce  a  high  capacity  flight  recorder 
suitable  for  scientific  experiments  aboard  shuttle.  There  are 
numerous  applications  for  a  2  gigabyte  recorder  so  initial  efforts 
have  concentrated  on  adapting  the  helical  scanning  technology. 

The  following  development  cycle  is  based  on  the  Exabyte  CTS. 

A.  Environmental  Specifications 

1 .  Vibration 


Mounting  the  tape  drive  requires  some  knowledge  of  the  shuttle's 
vibration  environment.  Vibration  and  induced  loads  will  vary 
depending  on  instrument's  mounting  location.  The  shuttle  was 
designed  to  carry  this  type  of  instrument  on  a  palette.  The 
mounting  rail  in  the  bottom  of  the  shuttle  bay  allows  slip 
movement  in  the  z  axis.  Palettes  can  therefore  flex  with  a 
natural  frequency  ranging  from  8  to  14  Hz.  To  prevent  scientific 
experiments  from  increasing  the  gain  of  a  palette,  the  natural 
frequency  of  an  experiment  must  be  no  less  than  35  Hz. 

The  following  table  illustrates  that  the  most  severe  induced  loads 
will  be  felt  in  the  z  direction. 


Table  1 . 

Enveloped 

Accelerations  (g) 

Lift  Off 

x-axis 

y-axis 

z-axis 

+  5 . 80 

+  1 . 55 

+13.15 

-9.00 

-1.55 

-12.28 

Landing 

x-axis 

y-axis 

z-axis 

+  7.47 

+  2.98 

+  8.19 

-7.70 

-3.03 

-5.12 

The  most  concerning  element  in  table  1  is  the  magnitude  of  the 
+13.15g  data  point  in  the  z-axis.  The  temporal  component  will  be 
obtained  to  determine  how  much  of  this  acceleration  can  be 
attenuated  by  isolators. 

The  sine  surveys  in  figure  CTS-1  were  performed  at  AFGL  and  point 
out  several  construction  points  that  can  be  improved. 
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x-axis:  This  curve  is  the  most  difficult  to  identify.  Removal  of 
the  door  and  tape  carriage  should  decrease  gain.  Fortunately,  the 
resonant  frequency  is  well  above  the  projected  natural  frequency 
of  the  vibration  mount  and  will  be  far  into  isolation  by  this 
point . 

y-axis:  The  lowest  of  the  resonant  frequencies  is  observed  in  the 

pendulum  like  swing  of  the  case  and  carriage.  Stiffening  by 

replacing  the  hanging  bracket  with  a  horizontal  strut  is  planned. 

z-axis:  The  vertical  axis  has  the  most  complex  response.  The 
first  peak  is  due  to  overhang  by  the  top  cover.  The  second  peak 
is  the  carriage  and  eject  mechanism.  Small  responses  are  believed 
to  be  other  tape  transport  mechanisms  and  the  backplane. 
Substitution  of  the  carriage  and  elimination  of  the  eject  feature 
will  be  the  first  modification. 

Analysis  of  a  first  try  vibration  isolation  can  be  found  in  figute 
CTS-2.  This  configuration  consists  of  off  the  shelf  parts  and  has 

not  yet  been  optimized  for  a  35  Hz  natural  frequency.  The 

isolators  are  mounted  at  four  points  located  radially  from  the 
equipment  center  of  gravity.  This  eliminates  second  order  peaks 
in  gain  at  higher  frequencies.  Isolator  materials  include  a  high 
damping  (C/Cc=0.15)  rubber. 

Figure  CTS-3  depicts  a  ruggedized  enclosure  /  vibration  mount  for 
the  tape  subsystem. 

2 .  Atmosphere 

A  peculiarity  with  magnetic  tape  is  that  it  requires  humidity.  If 
the  tape  is  allowed  to  dry  out,  the  ferrous  oxide  'flakes'  off  its 
Mylar  backing.  Loss  of  recording  media  can  be  prevented  by 
enclosing  the  recorder  inside  a  hermetic  box  with  a  controlled 
amount  of  moisture.  This  enclosure  is  also  required  because  of 
outgasing  that  may  be  present  from  the  recorder's  components. 

(This  commercial  recorder  has  not  been  assembled  from  vacuum 
compatible  parts).  The  flight  recorder  will  therefore  consist  of 
two  boxes,  a  vibration  mounted  box  within  a  hermetic  enclosure. 

B.  Electrical  Specifications 


The  tape  drive  (or  any  other  type  of  drive)  will  communicate  to 
the  DPU  via  a  Bridge  Controller  Board  (BCB).  The  BCB  will  be 
mounted  next  to  the  tape  drive  within  the  hermetic  box  and 
together  will  be  considered  a  flight  recorder.  The  BCB  will 
translate  the  SCSI  bus  into  an  efficient  serial  format  for  use  by 
the  DPU. 

Adding  a  BCB  isolates  the  quirky  housekeeping  overhead  of  specific 
drives  from  the  DPU.  This  black  box  approach  unburdens  the  DPU 
and  enhances  the  reliability  of  the  primary  objective,  detection 
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of  shuttle  charging.  Other  benefits  include  reduced 

EMI  &  EMC  through  the  use  of  coaxial  cable,  and  easier  swapping 

between  drive  technologies  that  use  a  SCSI  bus. 
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Isolator  Transfer  Function 


Ga'o  Squared 


1 
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FIGURE  QR3-14 
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SPREE  FLIGHT  RECORDER 
ENCLOSURE 


HIGH  VOLTAGE  OPTOCOUPLER 


The  first  custom  high  voltage  diodes  were  received  from  the 
manufacturer,  CSdc,  in  May.  These  diodes  were  designed  with  small 
bonding  pads.  The  small  pad  was  to  allow  IR  light  to  penetrate 
the  silicon  chip  and  generate  hole-electron  pairs.  The  diodes  did 
not  perform  as  expected  due  to  the  highly  doped  epitaxial  layer  on 
the  diode.  Silicon  is  normally  transparent  to  IR  light,  however, 
heavy  doping  in  the  p-layer  made  the  silicon  opaque. 

Our  next  attempt  will  etch  back  the  epitaxial  p-layer  to  expose 
more  of  the  diode's  depletion  region.  This  region  between  the  p- 
and  n-layers  is  where  photoconductivity  takes  place.  A  sample 
drawing  of  the  new  configuration  appears  in  figure  HV-1. 

Fiber  optic  faceplates  have  been  used  succesfully  as  a  light 
transmitting  medium.  Dielectric  tests  are  planned  to  demonstate 
the  optocoupler ' s  isolation  capabilities. 
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ROTARY  DRIVE 


In  order  to  increase  angular  coverage  of  the  270  degree 
analyser,  it  will  be  rotated  around  an  appropriate  axis. 

Since  lor  shuttle  missions  there  is  no  spacecraft  spin  to 
perform  this  function,  a  mechanical  system  must  be  developed 
to  rotate  the  sensors.  Oscillatory  motion  is  preferable  to 
continuous  rotation  because  it  eliminates  the  need  for  a 
rotary  interface  and  the  angular  coverage  may  be  adjusted. 
However,  oscillatory  m.otion  makes  much  greater  demands  on  the 
motor  since  it  requires  periodic  reversal  of  direction  of 
rotation . 

A  major  constraint  on  the  design  of  the  drive  system  is 
the  requirement  that  magnetic  fields  generated  by  the  drive 
motor  do  not  significantly  alter  electron  trajectories  inside 
the  analyzers.  Both  the  choice  and  location  of  the  motor  are 
important  in  this  respect. 

Several  approaches  to  this  problem  have  been 
investigated.  The  first  was  to  procure  an  integral  stepper 
motor,  speed  reducer,  and  rotary  platform  which  would  both 
support  and  drive  an  analyzer.  Such  a  device  is  sold  by 
Schaeffer  Magnetics,  who  call  it  a  "rotary  actuator."  The 
speed  reducer  utilized  is  an  especially  compact  type  known  as 
an  Harmonic  Drive.  From  discussions  with  Schaeffer,  it  was 
clear  that  this  approach  was  workable  and  relatively 
straightforward.  However,  the  cost  of  this  approach  was  a 
strong  deterrent  and  warranted  investigation  of  alternatives. 

Another  approach  considered  is  to  use  a  rotary  table 
mounted  on  a  large  diameter  ball  bearing  and  driven  by  Kapton 
belts  and  pulleys  from  a  motor/speed  reducer.  This  would 
have  the  advantage  that  the  drive  would  be  hollow,  ie,  a 
large  hole  would  exist  along  the  drive  axis  through  which  the 
cabling  could  be  routed,  minimizing  cable  length  and 
deflection.  Another  advantage  is  that  the  motor  could  be 
moved  as  far  away  from  the  sensors  as  necessary  in  order  to 
eliminate  magnetic  interference.  Sprocketless  Kapton  belt 
technology  is  well  proven  for  Shuttle  use  and  has  the  useful 
characteristics  of  vibration  isolation  and  zero  backlash. 

In  searching  for  appropriate  speed  reducers,  we 
discovered  that  compact  Harmonic  Drive  reducers  are  available 
from  their  manufacturer.  Harmonic  Drive  in  Wakefield,  Ma.,  at 
a  small  fraction  of  the  cost  of  the  Schaeffer  units.  One  of 
their  products  is  made  to  be  close-coupled  with  a  motor,  and 
has  a  planar  mounting  surface  at  its  output,  on  which  the 
sensor  may  be  mounted.  This  motor  -  reducer  combination  is 
equivalent  to  the  Schaeffer  unit,  with  the  added  flexibility 
that  the  motor  can  be  chosen  for  magnetic  characteristics  and 
can  be  mounted  further  away  from  the  sensor  if  necessary. 

This  system  is  reasonable  in  cost  and  relatively  simple  to 
implement  and  is  the  one  we  have  chosen  to  develop  and  test  a 
prototype  drive. 

Consideration  was  given  to  the  need  for  a  closed  loop 
control  system  using  an  optical  encoder  for  absolute  po' ition 
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feedback  to  the  motor  controller.  This  would  minimize  the 
possibility  of  positioning  errors  due  to  accumulated  step 
errors.  The  decision  was  made  to  build  and  test  the  prototype 
as  an  open  loop  system  with  optical  position  reference  points 
at  the  ends  of  the  scan  range.  This  approach  is  simpler  and 
the  optical  encoder  may  be  added  later  if  it  proves 
necessary . 
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Instrument  Simulation  Project 


2  -  D  Codes;  A  Personal  Computer  version  of  the  2-D 
codes  HEMIV  and  HEMIY  were  developed  to  give  us  simulation 
capabilities  independent  of  the  mainframe  computers  at  AFGL. 
This  was  done  using  the  Microsoft  Fortran??  V3.20  compiler 
with  mixed  success.  The  HEMIV  code  (using  elliptical 
trajectories  with  no  fringe  fields)  converted  easily  and 
produced  identical  results  on  286  and  386  level  PC's  with  or 
without  a  math  coprocessor.  HEMIY  (with  entrance  aperture 
fringe  fields  and  trajectory  tracing)  did  not  do  so  well. 

All  its  graphics  and  listing  capabilities  had  to  be  stripped 
out  and  the  grid  size  it  could  deal  with  minimized.  It  was 
then  divided  into  two  separate  programs,  the  first  half  to 
get  the  field  solution  and  the  second  to  do  the  geometric 
factor  integration.  The  development  was  done  on  a  286  class 
PC-AT  with  no  math  coprocessor  and  both  halves  ran  and 
produced  correct  results.  However,  attempts  to  run  it  on  a 
PC  with  a  math  co-processor  were  unsuccessful.  It  was 
tested  on  a  286  level  PC  and  on  a  386  level  PC  each  with  a 
8028?  math  coprocessor  and  it  bombed  in  both  cases.  The 
symptoms  were  different  on  the  two  machines  but  a  reasonable 
guess  would  be  that  the  problem  is  with  the  Microsoft 
interface  to  the  coprocessor  trig  routines.  A  decision  was 
made  to  table  the  project  until  a  better  Fortran  compiler 
can  be  acquired. 

3-D  Codes:  The  3-D  development  on  CLPH3D  continues. 
At  this  time  we  have  produced  a  field  solution  and  generated 
a  set  of  test  trajectories  and  traced  them  through  the  3-D 
grid.  A  comparison  between  the  central  plane  potentials  and 
those  from  the  equivalent  HEMIY  potential  solution  was  made. 
This  comparison  was  excellent  and  confirmed  that  the 
assumptions  made  for  the  2-D  code  fringe  fields  were 
correct.  We  are  now  evaluating  the  trajectories.  We  still 
have  not  implemented  the  interface  between  the  3-D  grid  and 
the  2-D  elliptical  code  through  the  interior  of  the 
instrument.  Also,  graphics  must  be  added,  and  the  geometric 
factor  integration  code  added.  The  grid  will  then  be 
expanded  in  the  third  dimension  and  the  CLPH3D's  simulation 
capabilities  explored. 
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PARTICLE  CORRELATOR 


Dr.  Paul  Gough 

12,  Briarcroft  Road 
Woodingdean 
Brighton,  BN2  6LL 
U.K. 


In  this  quarter  year  considerable  time  has  been  spent  looking  at  the 
problem  of  applying  particle  correlators  to  high  resolution 
spectrometers  of  the  type  that  might  fly  on  future  Shuttle  missions. 

One  of  the  greatest  problems  is  found  to  be  that  of  processing 
several  (~10)  fast  pulse  streams  simultaneously  without  having  massive 
duplication  of  electronics.  This  is  particularly  relevant  to  the  high 
frequency  (Mega-Hertz)  Buncher  mode  (electrons  only).  In  this  mode  the 
time  between  particle  arrivals  is  measured  by  hardware  in  units  of  a 
20  MHz  clock  in  each  of  the  ten  channels  separately.  It  is  not 
possible  to  perform  this  measurement  via  software  using  CPU  interrupts 
since  the  response  to  interrupts  is  of  the  order  of  a  microsecond, 
even  for  the  fastest  processors. 

The  hardware  is  required  to  measure  the  time  between  pulse  arrivals 
by  counting  the  20MHz  clock  and  to  present  those  measurements  to  the 
CPU.  This  initiates  a  software  increment  of  the  appropriate  memory 
location  to  generate  a  histogram  of  the  number  of  events  as  a  function 
of  time  separations  (64),  angle  zones  (10),  and  energy  levels  (32).  It 
is  assumed  that  these  histograms  are  read  out  via  telemetry  faster 
than  the  spectrometers  are  rotated.  (So  there  is  no  need  to  bin 
on-board  as  a  function  of  rotation  angle) . 

Each  pulse  stream  requires  one  six/seven  bit  counter,  a  six  bit 
latch,  a  J-K  flip  flop,  two  monostables,  and  four  Hand  gates.  The  J-K 
is  used  to  catch  overflow  (measurements  >64  units  of  50nS)  and  in 
those  cases  inhibit  the  incrementing  of  the  histogram.  With  this 
number  of  gates  per  pulse  stream  it  becomes  possible  to  consider  Gate 
arrays.  In  the  first  instance  it  was  thought  best  that  all  of  the 
Buncher  electronics  for  the  10  streams  could  fit  onto  one  Gate  Array. 
However  as  the  response  of  the  CPU  in  making  the  histogram 
incrementation  can  be  slow,  ten  separate  gate  array  units  might  be 
configured,  each  to  include  the  hardware  for  one  stream,  64  X  8  bits 
of  memory,  and  the  necessary  incrementation  logic.  Then  the  CPU  could 
read  out  the  full  Histogram  from  each  unit,  say  once  per  energy  level. 
In  this  way  the  detection  effeciency  (%  of  time  accepting  pulses  for 
measurement)  is  maximised.  Further  work  on  this  aspect  are  awaiting 
discussions  in  june  with  Amptek  who  have  more  experience  with  gate 
array  design. 
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of  a  few  tens  of  mS  (or  say  the  energy  step  dwell  time)  into  four 
equal  time  bins  corresponding  to  wave  phases  0-90,  90-180,  180-270,  & 
270-360  degrees.  If  there  is  no  significant  modulation  then  the 
variation  of  these  four  sums  should  be  of  the  order  of  that  expected 
by  Poisson  statistics,  (i.e.  differences  of  the  order  of  the  square 
root  of  the  average  total),  it  is  a  simple  calculation  to  find  the 
difference  between  the  greatest  and  the  smallest  and  divide  that  by 
the  square  root  of  the  average.  This  gives  a  direct  measure  of  the 
significance  of  any  modulation  present  in  the  particles  at  the  beam 
modulation  frequency.  Data  can  be  very  compressed  with  the  CPU  only 
transmitting  significant  data  (  >  3  sigma).  Furthermore  the  CPU  can 
use  the  knowledge  of  where  in  velocity  space  the  modulation 
significance  is  highest  to  concentrate  the  resolving  power  to  that 
region . 


WORK  PLAN 


The  scheduled  research  falls  into  a  number  of  specific  topics  listed 
here : 

1)  An  early  design  step  is  to  choose  the  number  of  resolution  bins  in 
frequency,  energy,  zone  angle,  and  rotation  angle  for  each  of  the  two 
particle  species.  For  a  given  data  rate  this  naturally  implies  a 
tradeoff  between  these  measured  parameter  resolutions  against  temporal 
resolution.  It  also  sets  the  required  m.emory  size  for  holding  these 
average  ACFs . 

2)  Study  CPU  types  available  and  choose  an  appropriate  processor.  At 
this  time  the  transputer  looks  like  being  a  likely  choice.  In  the  near 
future  an  auroral  rocket  will  use  a  group  of  4  transputers  to  replace 
36  microcomputers.  The  transputer  has  the  advantage  of  true  parallel 
processing  and  can  be  arranged  in  fault-  tolerant  configurations. 

3)  Verify  that  the  chosen  CPU  type  can  provide  the  processing  required 
via  a  software  simulation.  In  the  case  of  transputers  this  simulation 
will  determine  the  number  of  transputers  required. 

4)  Design  the  hardware  block  components.  Particular  problem  areas  are: 
a)  Fast  pulse  timing  for  the  Buncher  mode  on  ten  data  streams,  (b) 

Fast  memory  accessing  by  transputers,  (c)  Effective  division  of  the 
modulation  reference  signal  into  four  equal  time  bins. 

5)  Liase  with  AMPTEK  over  the  specific  hardware  design  with  a  view  to 
AMPTEK  constructing  the  particle  correlator  cards. 

6)  Liase  with  A.MPTEK  over  the  flight  software.  It  is  expected  that  the 
software  will  be  a  joint  production. 

7)  Design  a  checkout  system  that  can  effectively  test  the  required 
flight  functions. 

8)  Support  AMPTEK  in  constructing,  testing,  checkout,  launch 
preparations,  in  flight  monitoring,  and  subsequent  post  flight  data 
interpretation . 
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EXPECTED  FREQUENCIES  &  TECHNIQUES 


As  stated  above  the  type  of  particle  correlation  method  depends  on 
the  expected  frecjuencies  of  modulation: 

(a)  Natural  1-10  MHz  modulations 


Electrons  (only)  are  expected  to  undergo  interactions  with  waves  at 
these  high  frequencies:  the  local  upper  hybrid  frequency,  Fuh  (1-7 
MHz)  and  local  electron  gryof requency ,  Fg  (1.4MHz-  orbit  location 
dependant)  and  its  low  harmonics.  Bernstein  or  Electron  Cyclotron 
Harmonic  Waves  are  to  be  expected  at  (n  +  1/2)  Fg  and  at  Fuh. 

At  these  high  modulation  frequencies,  approximately  ten  times  higher 
than  the  highest  expected  count  rates,  the  one  bit  'buncher'  mode  is 
the  ideal  merhod  of  correlation.  This  method  involves  measuring  the 
time  of  arrival  between  successive  electrons  in  units  of  say  a  20MHz 
clock.  A  histogram  of  the  number  of  measured  separations  as  a  function 
of  separation  units  (50nS)  is  accumulated  over  periods  typically  of 
lOS  or  more  depending  on  the  average  count  rate.  This  method  is 
directly  equivalent  to  averaging  many  one  bit  correlation  functions, 
each  with  only  two  bits  set.  It  is  much  easier  to  produce  in  hardware 
and  there  is  neglible  loss  in  detection  efficiency  since  the 
probability  of  obtaining  more  than  two  bits  set  in  say  a  one  bit 
correlator  data  take  of  64  samples  at  20MHz  is  neglible  for  typical 
count  rates  of  a  few  times  100,000  per  second. 

This  technique  has  been  used  successfully  to  identify  modulations  at 
mega-  Hertz  frequencies  present  on  the  natural  auroral  beams  above 
auroral  arcs. 

(b)Natural  O-lOKHz  modulations 


This  lower  frequency  range  contains  both  the  lower  hybrid  frequency 
(1-lOKHz)  and  the  ion  gyrof requencies  (  <  IKHz).  Naturally  occurring 
modulations  at  these  frequencies  are  expected  both  in  ions  and 
electrons . 

Average  Auto-Correlation  Functions  (ACFs)  can  be  generated  in 
software  from  say  groups  of  128  sant'le  takes  at  20KHz.  On  AMPTE  UKS 
this  technique  was  further  simplified  by  converting  each  set  of 
samples  first  to  one  bit  by  comparing  each  sample  with  the  average  of 
that  set.  One  bit  ACFs  were  rapidly  calculated  and  averaged.  At 
kilo-Hetrz  frequencies  a  NSC800  processor  can  achieve  >90%  detection 
efficiencies  (%  of  time  taking  samples).  The  disadvantage  of  this 
technique  is  that  it  is  only  sensitive  to  the  strongest  modulation 
present.  With  the  higher  computing  powers  now  available  (e.g. 
transputers)  full  multi-bit  ACFs  should  be  possible.  There  are  also 
dedicated  ACF  chips  but  microprocessor  based  designs  offer  greater 
flexibility  given  the  range  of  instrument  modes  and  frequency-  energy- 
2-D  angular  binning  needed  for  the  results. 

(c)  Artificial  Beam  Modulation 


As  mentioned  above  the  active  electron  and  ion  beams  might  be 
deliberately  modulated  at  kilo-Hertz  frequencies.  If  the  frequency  and 
phase  of  this  reference  signal  is  present  on  board  then  the 
spectrometer  detected  particle  samples  can  be  directly  cross 
correlated  to  look  for  this  modulation.  In  principle  a  cross 
correlator  could  be  achieved  as  in  (b)  above  but  a  specialised 
technique  has  many  advantages. 

If  the  reference  signals  can  be  generated  at  four  times  the 
modulation  frequency  by  PLL  then  counts  can  be  summed  up  over  periods 
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When  particles  are  in  velocity  resonance  with  waves,  either  particle 
energy  is  transferred  to  the  waves  resulting  in  wave  amplitude  growth 
and  particle  decceleration ,  or  wave  energy  is  transferred  to  the 
particles  resulting  in  particle  acceleration  and  wave  damping.  In 
either  case  particles  in  velocity  resonance  will  become  phase  bunched 
with  the  wave.  This  is  observable  as  fluxes  of  particles  at  the 
resonant  velocity  (or  appropriate  energy  level)  being  modulated  in 
time  at  the  wave  frequency.  Measurements  of  these  particle  modulations 
as  a  function  of  frequency,  energy  level,  and  particle  direction  of 
motion,  allow  us  to  identify  directly  by  experiment  those  regions  of 
particle  velocity  space  contributing  to  wave  growth  and  wave  damping. 

Active  beam  experiments  with  beams  actively  modulated  at  kilo-Hetrz 
frequencies  allow  for  the  beam  electrons  or  ions  to  be  tagged  or  to 
stimulate  natural  emissions  at  these  frequencies.  .Any  beam  electrons 
returning  to  the  Shuttle  bay  should  also  be  modulated  at  this 
frequency  and  identified  directly  by  a  particle  correlator. 

Furthermore  any  beam  electrons  which  have  undergone  acceleration  or 
decceleration,  or  natural  electrons  which  have  interacted  with  them 
via  wave-  particle  interactions,  will  also  be  modulated  at  the  beam 
modulation  frequency  and  identified  as  such. 

There  are  several  methods  which  can  be  employed  by  which  particle 
modulations  can  be  measured.  The  choice  of  method  depending  mainly  on 
the  frequency  range  of  interest  and  to  a  lesser  extent  on  other 
information  available  on-board  at  the  particle  correlators  (i.e 
mission  specific).  For  a  realistic  design  a  specific  particle 
spectrometer  is  assumed  here. 


EXAMPLE  APPLICATION 


For  initial  hardware  design  purposes  it  is  assumed  that  there  are 
two  pre-existing  particle  instruments,  each  a  combined  electron  and 
ion  spectrometer.  It  is  further  assumed  that  these  instruments  are  of 
the  270  degree  type  of  electrostatic  analyser  with  each  instrument 
having  a  field  of  view  100  x  10  degrees  (  in  the  form  of  ten  detection 
zones  each  of  10  x  10  degrees  spread  around  in  a  fan).  It  is  also 
expected  that  these  two  instruments  are  mounted  in  such  a  way  that 
thay  can  be  rotated  about  an  axis  perpendicular  to  the  plane  of  these 
fan  field  of  views.  A  resulting  total  field  of  view  is  achieved 
approaching  the  maximum  2pi  available  to  a  Shuttle  bay  mounted 
instrument  (i.e.  subject  to  limits  imposed  by  shaddowing  from  other 
instruments  in  the  bay). 

Active  electron  and  ion  beam  experiments  usually  allow  for  the 
frequency  modulation  of  the  beam  at  kilo-Hertz  frequencies.  It  is 
assumed  in  this  design  study  that  a  signal  corresponding  in  frequency 
and  phase  to  this  beam  modulation  frequency  is  available  on-board  for 
use  as  a  reference  signal  by  the  particle  correlator. 

As  there  are  two  identical  particle  instruments  (each  electrons  and 
ions)  it  is  planned  to  have  identical  correlator  units  in  each  one. 
Thus  the  inputs  available  to  each  of  these  correlator  units  are  20 
pulse  streams  (  2  species  x  10  zones  ) ,  one  beam  modulation  reference 
signal,  and  some  instantaneous  monitor  of  the  spectrometer  energy 
level. 
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SENSOR  DESIGN 


Hemispheres  necessary  for  the  building  ot  a  nested  tri-quadraspheric 
electrostatic  analyzer  (Q^ESA)  were  ordered  and  received  this  quarter.  These 
hemispheres  were  manufactured  to  our  specified  outer  and  inner  diameters  from 
type  6061-T6  aluminum.  The  hemispheres  have  been  stress  relieved  as  part  of  their 
fabrication.  They  are  specified  to  be  within  r.OOl  inch  inner  diameter  and  +.000/- 
.005  of  their  outer  diameter.  Their  wall  thickness  is  .050. 

This  is  material  acceptable  for  spaceflight  hardware.  The  stress  relief  and 
tolerances  are  necessary  to  insure  that  the  machined  and  mounted  quadraspheres 
are  mutually  concentric.  A  high  degree  of  concentricity  minimizes  fringe  fields  that 
may  distort  an  Q^ESA’s  responses. 

A  vacuum  chamber  thermal  control  shroud  was  installed  and  tested  in  the 
vacuum  chamber.  At  first  the  shroud  performed  properly,  but  after  several  cold/heat 
cycles  it  developed  a  leak.  This  leak  must  be  located  and  repaired  before  the  system 
is  used  to  test  any  contamination  sensitive  instruments  such  as  the  Q^ESA.  The 
thermal  vacuum  chamber  testing  is  necessary  to  verify  new  designs  to  insure  they  can 
survive  a  space  environment. 

A  visit  was  made  to  Martin  Marietta  in  Denver  Colorado  to  discuss 
mechanical,  electrical,  and  thermal  interface  problems  typical  to  a  space  shuttle 
flown  payload.  Martin  is  an  integration  contractor  for  several  shuttle  payloads  and 
has  extensive  experience  in  developing  design  criteria  for  successful  shuttle 
experiments.  Specifics  of  the  environment,  including  vibration  levels,  operating 
thermal  ranges,  typical  electrical  interfaces,  and  mandatory  safety  requirements  were 
covered.  Development  of  a  prototype  space  qualified  Q^^SA  requires  not  just 
meeting  scientific  performance  criteria,  but  in  also  meeting  environmental 
specifications  of  spaceflight. 

Amongst  many  useful  subjects  covered  in  the  Denver  meeting  was  the 
possible  use  of  a  "cold  plate"  to  provide  thermal  control  to  shuttle  bay  experiments. 
This  device  attaches  to  a  controlled  temperature  fluid  circulating  system  in  the 
shuttle  which  regulates  the  temperature  of  the  cold  plate  to  a  +2®C  to  +50°C 
temperature  range.  Assuming  the  use  of  such  a  device  simplifies  the  development 
and  design  of  the  sensor  hardware  since  it  greatly  reduces  the  possible  thermal 
extremes  that  an  instrument  may  see  in  spaceflight. 

Th“  area  of  a  "cold  plate"  is  limited.  The  mounting  feet  of  the  tape 

recorder  were  reoriented  to  accommodate  mounting  to  the  smaller  surface  area. 

(See  the  attached  drawing.) 
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C-0094  Date:  07-22-88  bwg.  No.  D:\SPREE\RTMD 


DATA  PROCESSING  UNIT  FOR  270  DEGREE  NESTED  SENSORS 


The  parts  selection  and  design  for  the  CPU  boards  is  now  complete.  A 
prototype  CPU  board  was  constructed  to  verify  the  proper  operation  of  the  memory 
chips  and  several  other  circuits.  The  PC  board  layout  is  now  underway  and  the  first 
PC  boards  should  be  available  for  testing  during  the  next  quarter. 


The  following  parts  are  being  used  for  the  CPU  board: 
Part  # _ Description 


80C86-2 

82C84A 

82C59A 

82C54 


8  MHz  CMOS  microprocessor 
Clock  generator 
Interrupt  controller 
Timer/counter 


HM-6j642 

HM-6617B 

EP320 

74AHCT373 

74HC245 

74FCT244 

IDT7 130/7 140 


8k  X  8  static  RAM 
2k  X  8  CMOS  PROM 
Altera  EPLD 

8-bit  latch 
8-bit  transceiver 
8-bit  buffer 

Master/slave  dual-port  RAM 


A  preliminary  PC  board  layout  was  performed  for  the  D400  accumulator 
board  to  verify  that  the  proposed  board  size  of  9.75  '  x  6.25"  is  appropriate.  It  is,  and 
the  final  D40d  board  layout  will  be  completed  during  the  next  quarter. 


Included  is  a  partial  block  diagram  of  SPACE  (Spacecraft  Particle  Correlator 
Experiment).  Only  half  of  SPACE  is  shown  in  the  block  diagram,  the  other  half  being 
identical.  The  operation  of  the  various  modules  has  been  described  in  previous 
quarterly  reports. 

Not  shown  in  the  block  diagram  is  the  control  CPU.  SPACE  will  be  controlled 
by  an  80C86  with  a  full  complement  of  suppon  chips.  Communications  between 
SPREE  and  SPACE  will  be  performed  through  a  pair  of  dual-port  RAMs. 

The  microcontroller  modules  for  SPACE  will  consist  of  an  80C31 
microcontroller,  an  HM6516  static  RAM  and  an  87C64  EPROM.  Not  shown  on  the 
block  diagram  is  a  redundant  microcontroller  module,  which  will  automatically  switch 
on-line  in  case  one  of  the  primary  microcontroller  modules  fails. 


The  high-frequency  buncher  modules  will  contain  an  Altera  EP1800  EPLD, 
an  HM6516  and  other  support  logic.  This  module  has  been  prototyped  with  discrete 
components  and  is  now  being  fit  into  programmable  logic  devices  to  reduce  board 
space  requirements. 
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SPREE  FLIGHT  RECORDER 


Amptek  has  decided  to  use  helical  scanning  tape  technology  for  the  flight 
recorder.  A  vendor.  Exabyte  of  Boulder  CO,  has  been  chosen  to  supply  commercial 
grade  hardware.  Amptek  will  repackage  the  drive  for  space  flight. 

Two  members  from  Amptek  visited  Exabyte  in  Au^st  to  discuss  the 
recorder’s  intended  application.  After  talking  to  several  of  the  recorder’s  designers, 
we  are  confident  that  this  is  a  wonhy  device.  Information  was  also  obtained  for 
methods  to  enhance  the  drive’s  reliability. 

The  flight  recorder,  currently  designated  FR20.  will  store  over  one  billion 
bytes  of  data  during  the  SPREE  mission.  For  redundancy,  two  recorders  will  be 
flown  side-by-side. 

Present  work  involves  the  design  of  a  hermetic  box  for  the  drive  mechanism. 
Magnetic  tapes  requires  humidity  to  retain  their  magnetic  coating  so  the  enclosure 
will  contain  a  controlled  environment  The  planned  Bridge  Controller,  BC,  will 
modify  the  internal  air  temperature  of  the  box  before  the  tape  mechanism  is 
initialized. 

Signet  Tool  &  Engineering  has  been  contacted  to  manufacture  the  enclosure. 
We  are  planning  to  maclune  the  base  from  a  single  block  of  aluminum.  The  cover 
will  also  be  made  from  a  single  sheet-  All  external  bolts  will  be  captured. 
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HIGH  VOLTAGE  OPTOCOUPLER 

The  high  voltage  photodiode  has  evolved  a  new  design.  The  previous 
approach  with  single  junctions  has  been  abandoned.  A  more  compact  and  efficient 
diode  stack  is  now  being  considered.  Delivery  of  sample  parts  is  expected  at  the  end 
of  September. 

A  May  plant  trip  to  CSdc  in  New  Jersey  helped  advance  the  optocoupler’s 
design.  Observation  of  manufacturing  techniques  and  capabilities  suggested  that  a 
stacked  junction  diode  would  be  a  much  better  component. 

A  wire  frame  model  of  the  proposed  optocoupler  follows.  Two  prominent 
features  are  the  high  voltage  outlets  and  the  method  for  heat  removal.  For  high 
voltage,  hollow  pins  are  cast  into  the  optocoupler’s  epoxy  body.  Their  location  is 
chosen  for  highest  isolation.  Connections  can  be  soldered  to  the  pins.  Heat  removal 
is  accomplished  by  two  tin  plated  copper  bars  on  either  side  of  the  body.  The  bars 
are  soldered  to  heat  sinks  on  a  printed  circuit  board.  Leads  for  the  LEDs  are  also 
soldered  to  the  p.c.  board. 

The  first  prototype  will  have  a  blocking  voltage  of  12kV  and  couple  80uA  of 
output  current  for  20mA  of  input  current.  The  unit  svill  measure  0.700"  x  0.250"  x 
0.200" 
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ISOMETRIC 


HIGH  VOLTAGE  OPTOCOUPLER 

FIGURE  QR4-3 


SPACECRAFT  PARTICLE  CORRELATOR  EXPERIMENT 


Work  of  consultant:  Dr.  Paul  Gough.  Brighton.  UK 


An  early  design  step  has  been  to  choose  the  number  of  resolution  bins  in 
frequency,  energy,  zone  angle,  and  rotation  angle  for  each  of  the  two  particle  species. 
This  sets  the  required  memory  size  for  holding  these  average  ACFs  and  the 
telemetry  requirement.  At  the  same  time  Electron  Gun  Pulse  patterns  have  been 
specified.  In  particular  H.F.  buncher  and  pulse  patterns  have  been  designed  to 
enable  an  "Electron  Radar  '  mode  of  operation. 

Considerable  time  was  spent  studying  CPU  types  available  with  eventual 
choice  of  iiPAK607C31  as  an  appropriate  processor  that  combines  processing  power 
with  large  memory  (PROM  +  RAM)  in  one  40  pin  DIL. 

Over  this  quarter  the  main  task  has  been  to  design  the  hardware  of  the  fast 
pulse  timing  for  the  Buncher  mode.  I  have  constructed  and  fully  tested  a  buncher 
breadboard  unit  formed  out  of  74HC  parts  that  is  fully  microprocessor  bus 
compatible.  This  design  is  now  being  studied  by  AMPTEK  with  a  view  to  reducing  a 
large  part  of  the  circuit  onto  a  gate  array  chip.  A  selectable  clock  will  be  included  to 
give  a  variety  of  frequency  ranges,  equivalent  to  distance  ranges  in  the  "Electron 
Radar"  mode. 

This  work  is  part  of  the  overall  liaison  with  AMPTEK  over  the  complete 
hardware  design  with  a  view  to  AMPTEK  constructing  the  particle  correlator  cards.. 
The  software  of  the  chosen  microprocessor  type  is  now  being  studied  with  a  view  to 
the  software  writing  being  a  joint  effort  with  AMPTEK. 
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JX.  jBi  bnrmi. 


INSTRUMENT  SIMULATION  PROJECT 


An  advanced  Fortran  compiler  for  Personal  Computers  was  purchased  to  run 
the  HEMIV  and  HEMIY  2-D  simulation  codes.  This  compiler  had  to  be  fully  main 
frame  compatible  to  accommodate  the  developed  and  proven  simulation  codes.  A 
Ryan-.McFarland  Fortran  was  selected.  This  product  has  been  implemented  and 
HEMIV  and  HEMIY  codes  are  now  executing  on  our  PC.  Attached  are  listings  of 
the  current  state  of  these  two  programs. 

Instructions  for  building  a  input  data  file  for  HEMIY  were  written.  This 
allows  a  designer  to  model  and  test  a  design  using  the  HEMIY  code  without 
involving  a  programmer  to  execute  the  program.  The  development  of  this  sort  of  in- 
house  "user  friendly"  code  was  a  major  objective  of  our  simulation  development 
program.  A  copy  of  these  instructions  is  enclosed. 

3-D  code  development  on  CLPH3D  has  continued.  This  code  still  requires  a 
main  frame  processor.  Output  results  between  it  and  the  PC  run  HEMIY  are  being 
made  and  to  an  extent  they  are  being  used  to  cross  prove  each  other.  The  fringe  field 
solutions  of  CLPH3D  are  very  accurate  and  can  be  used  to  test  2-D  approximations 
used  on  HEMIY.  On  the  other  hand,  the  trajectory  analysis  done  by  HEMIY  is 
being  used  to  confirm  the  output  of  the  more  computational  rigorous  CLPH3D.  The 
interface  between  the  3-D  grid  and  2-D  elliptical  code  has  been  started.  A  graphics 
display  and  output  is  being  developed  to  help  interpret  and  understand  the  results. 
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HEMIY  DATA  FILS  DOCUMENTATION 


The  first  step  in  using  the  sinaulation  code  developed  by 
AMPTSK  is  to  construct  a  data  file  defining  a  2-D  grid,  the 
instrument  walls.  a  set  of  apertures,  and  the  location  of  all 
relevant  fixed  voltages.  The  data  file  used  by  HEMIY  specifies 
the  above  information  in  an  abbreviated  format.  In  order  to 
define  it,  first  consider  Figure  1.  Figure  1.  illustrates  a  12  x 
8  grid  which  is  of  adequate  size  to  illustrate  all  the  relevant 
features  modelled  by  HEMIY.  The  two  plates  of  the  instrument, 
labeled  RI  and  RO  are  shown  in  dark  black,  the  four  apertures, 
labeled  A1  through  A4  are  shown  in  red,  and  the  remaining  fixed 
voltage  line  segments,  FI  through  F6  are  shown  in  blue. 

The  data  file  corresponding  to  Crude  Instrument  consists  of 
fifteen  lines  of  information  as  shown  in  Table  1.  below: 


Cl: 

DATA  FILE  FOR  CRUDE 

INSTRUMENT. 

AUGUST  1988 

C2; 

12  8  2 

8 

C3:4 

,4,4,5  6, 6, 3, 6 

10, 

10,3,6 

12, 

12,3,6 

C4:6 

1,2.0  2,5.5 

3,6.0 

6,8 

.0  7,8.5  8,11 

C5;3 

1,90.0  6,90.0 

12, 

-180.0 

C6:3 

l,-5.0  6,0.0 

12,0.0 

C7:3 

l,-5.0  6,0.0 

12,0.0 

C8: 

150.0  6,12,3,3 

4 

C9:- 

137.0  6,12,6,6 

5 

CIO: 

0.0  4,12,7,7 

6 

Cll: 

0.0  4, 4, 1,4 

6 

C12: 

0.0  4, 4, 5,7 

6 

C13: 

0.0  12,12,1,3 

6 

C14: 

0.0  12,12,6,7 

6 

C15: 

-98.0  7, 7, 3, 6 

9 

TABLE  1. 


Card  labels  'Cn:'  have  been  added  for  reference  and  are  not  part 
of  the  real  data  file.  Cl  is  an  80  column  header  card;  C2  gives 
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vhe  four  values  MX,  NY.  NL ,  AND  NV ;  C3  gives  the  four  apertures; 
C4  defines  the  vector  quantity  R';  C5  defines  'AL?';  C6  defines 
' X3 ' ,  and  C3  through  C15  specify  fixed  voltage  line  segments. 

HEMIY  constructs  an  NX  x  NY  distorted  rectangular  grid  with 
NL  logical  levels  Each  grid  point  has  (x,y)  coordinates  based  on 
the  three  vector  quantities  R,  ALP,  and  X3.  Note  the  shorthand 
used  in  C4.  C5  and  C6.  For  example  the  R  vector  is  specified  by 
six  number  pairs  each  consisting  of  an  index  and  an  R  coordinate. 
HEMIY  calculates  R(4)  and  R(5)  by  linear  interpolation. 
Similarly,  only  three  pairs  are  needed  to  define  ALP  with  ALP(l) 
through  ALP(6)  having  the  value  90.0  and  ALP(7)  through  ALP(ll) 
containing  values  at  45  deg.  intervals  from  45.0  through  -135.0. 
X3  varies  linearly  from  -5.0  at  XB(1)  to  0.0  at  XB(6)  and  is  zero 
for  X3(7)  through  XB(12).  The  grid  coordinate  algorithm  should  be 
fairly  obvious  by  now.  HEMIY  treats  R  as  a  Y-coord  and  XB  as  a  X- 
coord  whenever  ALP  is  90.0.  It  treats  R  and  ALP  as  a  polar  grid 
whenever  ALP  is  not  90.0. 

You  may  have  noted  that  nothing  has  been  said  about  C7  and 
about  what  NL  does.  That  is  because  at  this  time  they  do  nothing 
but  must  be  included  in  the  data  file  for  historical  reasons. 
When  the  data  file  structure  was  originally  designed,  I  had  in 
mind  modelling  the  box  containing  the  instrument.  This  would  have 
required  a  non=poiar  region  between  RO  and  the  box.  Since  this 
was  just  a  cosmetic  enhancement,  it  was  never  implemented.  Note 
that  the  value  of  NL  is  ignored,  but  that  C7  must  be  identical  to 
C6  . 

C3  defines  the  four  apertures  shown  in  Figure  1.  Each 
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aperture  is  defined  by  a  set  of  four  integers,  IXl , 1X2 , JYl , JY2 . 
Hence  the  aperture  is  the  line  segment  connecting  points  (1X1, JYl) 
and  (IX2,JY2).  A1  is  the  input  aperture,  and  IIEMIY  uses  the  A1 
coordinates  to  define  the  injection  plane  and  the  size  of  the 
aperture,  A2  defines  RI  and  RO,  A3  is  the  130  degree  aperture 
which  is  no  longer  used  and  A4  is  the  exit  aperture  which  is  also 
ignored  at  this  time. 

The  remaining  cards  are  the  voltage  cards.  There  are  NV  of 
them  and  the  first  two  are  the  inner  and  outer  plate  voltages. 
The  format  is  simple:  a  voltage  value,  four  integers  specifying  a 
line  segment  using  the  C3  notation  and  an  integer  flag.  They  were 
intended  to  be  order  independent  since  the  flag  is  sufficient  to 
define  them,  but  that  is  not  guaranteed.  The  flag  values  for  the 
plates  must  be  4  and  5  and  the  other  fixed  voltages  should  be 
flagged  with  a  6.  CIS  gives  the  interface  plane  between  the 
trajectory  code  region  and  the  region  where  the  elliptical 
trajectory  is  applied  and  it  must  be  flagged  with  a  9.  Note  that 
the  voltages  specified  on  C8 ,  C9 ,  and  CIS  are  replaced  with 
appropriate  values  by  HEMIY.  The  ground  voltages  on  CIO  through 
C14  are  used  as  supplied  and  hence  should  be  zero. 

One  final  note  on  the  data  file  format.  Cl  is  an  ASCII 
label  and  is  read  with  a  Fortran  (A80)  format.  C2  must  be  (4IS) 
since  there  are  some  other  diagnostic  parameters  read  from  C2  that 
should  all  default  to  zero.  The  remaining  card  images  are  read  as 
Fortran  free  form,  meaning  that  their  values  may  be  separated  by 
commas,  one  or  more  spaces,  or  a  carraige  return,  but  that  there 
must  be  the  correct  number  of  values. 
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oo 


'P.CGRAM  HEMIV 

IMPLICrT  REALMS  (A-rl.Q-Z^ 

d-IARACTSS  TEST(lOl.lOl)  ,rCHAR(10)  , -TMP  (  ' *  1 4 
DIMENSION  :T0T(  lOL.lOl )  ,.'ISUM(21  1  ,:StJM(  101  !  ,  JSUM(  101 )  , 

-RHCD( 15) ,LSUM(101) ,MSUM(101) ^IJEOTtZl) ,NACT;21) 

2IMENSICN  XX(756)  .ZY(7S5)  ,MCON(  7 , 1394  )  .W(7S5)  ,MBDX(2,130)  , 
-ANEL(1294)  ,3K(44,756)  ,CHD(  1394)  /■/?  ( 1  37  )  ,  ?P  '  756  )  ,RR  (  756  )  ,IAP(16)  , 
-:N(41 ,JCCN(26) ,AINV(3,5) ,rAU(1294) ,IBK(4) ,WW( 17.18) 

COMMON  /ORB/  XDOT,YDOT,XP//P,NTSAV,?IMIN,NOWEr,.NOWSO 
COMMON  /INST/  RIP. R0P.3IA.R0A. SNOT, £K 
OATA  ISUM.JSUM.OSUM.MSUM  /404»0/ 

DATA  ICHAR/ '  •,'2',‘3','4',*5’,'6','7','9','9','9'/ 

DATA  TAU  /1394*0./ 

DATA  TEST  /10201*’  '/ 

DATA  NEOT.NAOT  /42*0/ 

DO  10  1*1.41 
DO  10  J-1,41 
ITOT(I.J)-0 
10  CONTINUE 

DO  20  N-1,21 
NSUM(N)*0 
20  CONTINUE 
NTRJ-0 
PI*ACCS(-1.) 

RAD-PI/130. 

THETA-90. 

THETA-THETA »HAD 
ST-3IN( THETA) 

Miao-0 
C  MTHD-O 

REWIND 7 
RSWIND4 
RFWIND8 
REWINDll 

SET  QUADRATURE  LIMITS 
ME-41 
MP-41 
MR-11 
MA-11 

C  CALL  SSWTCH(6,ISW6) 

ISW6-1 

IFPLT-1 

I?(ISW6.EQ.l)  THEN 

PRINT*. 'M180,MB,MF,MR..MA  -  '  .MISO.ME.MP.MR.MA,  '  ENTER  NEW  VALUE 

vS' 

READ(5,*)  Mia0.ME,MP.MR,.MA 

PRINT*. 'M180,ME,MP,MR..MA  -  '  ,M180  .ME..MP.MR..MA 
IFPLT-2 

IP(M180.GT.3)  THEN 
M180-M180-4 
C  MTHD-1 

SNDIF 
ENDIF 
C 

C  DEFINE  GEOMETRY 
ENOT-1000. 

CALL  GECMB(VOUT.VINN.ai?,HCP,RIA.RCA.SNOT.NX,.NY,.NBDY,XX.YY.VV.CHD. 
-VP , MBDY . PP . RR . LAP . AEND . IPRT . IBK . NBK ) 

.'rrOT-NX*MY 
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on  oo  on  o  o  non  no 


.'fELa  NX-  :  ;  *  'NY-1 )  »2 

^ETINirrON  OF  ?EAT2  APERTURE  RADII 

RB=  ?.I?*RC?^/2 
ORF’-RCP-RI? 

OIA«RCA-RIA 
?A»0IA/; . 

EAl^OIA 


SET  130  OEG  APERTURE  ?.\RAMETSRS 
IF'MiaO.EQ.l)  THEN 
SLIT  OF  CONSTANT  DELTA  R 
RC-2AL/2. 

WRITS ( 6, ' (IQX, ■ 'APERTURE  AT  130  IS  SLIT. 
-  ) ' )  SAL 

ELSE  IF!MiaO.SC-2)  THEN 
CIRCLE  OF  RADIUS  RC 
RC»RA 


WIDTH  ' ,F8.3, ' ' 


WR*.S(S,  ( Al, 6 X, ’ 'APERTURE  AT  130  IS  CIRCLE.  RADIUS 
CM'')')  ■  '.RC 

ELSE 

DEFAULT  IS  NO  APERTURE 
M130-0 

RC«'RCP-RI?)/2. 

WRITE(S, ' (SX, • 'NO  APERTURE  AT  180  DEGREES'')') 

SNDIF  )  I 


m>  ' 


,P8 


DEFINITION  OF  U  QUADRATURE  PARAMETERS 
UIP-RB/RI? 

UCP-aa/RQP 
UIA-R3/RIA 
UOA-RB/RCA 
UINC»(UIA-UOA)/MR 
UST-UOA+OlNC/2 
IF(MR.LE.l)  THEN 
.MR»1 
UST-i. 

SNDIF 


DEFINITIONS  OF  ALPHA  QUADRATURE  PARAMETERS 
ALF31— ASIN(RA/RB  ) 

ALF32»-ALF31 
ALFINC- ( ALF32-ALFB1 ) /MA 
ALFAST-ALF31+ALFINC/2 
IF(MA.LE.l)  THEN 
MA-1 

ALFAST-0. 

SNDIF 

DEFINITION  OP  EK  QUADRATURE  PARAMETERS 
EXl-0.85 
EK2-1.15 

SKINC* ( EX2-2K1 ) /ME 

SKST»EKlt-EKINC/2 

DE-(VINN-70UT)/21000. 

SB-1 .-10.*DE 
IF'.ME.LS.l)  THEN 
ME-1 
E.Kl-1. 

SK2-1. 


CM'  ' 


3. 
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iXINC-l . 

SXST-1 . 

SNDir 

c 

C  DEFIMiriCN  OF  PHI  QUADRATURE  PARAMETERS 
PHi:=i-9.Q/ST 
PHIO^IS.O/ST 
PHI1*?HII*RAD 
?HI2»?KI2*RAD 
PKINC* ( PHI 2- PHI 1 ) /MP 
PHIST-PHIl+PHINC/2 
IF{MP.2E.l)  THEN 
MP-L 
PHIl-0. 

PHI2-0 
PHINC-1 ./ST 
PHIST-0. 

ENDI7 

ASUM-0 

RIA-RB-RC 

ROA-RB-t-RC 

UIA-RB/RIA 

UOA-RB/ROA 

C 

C  DEFINE  CONNECTIVITY 

CALL  CNEC{NX.NY,NEL,NBDY,NBDW.MCON,MBDY,ANEL,XX.yY,VV,VP) 

C  SOLVE  FOR  POTENTIALS 
DEB YE- 0. 

CALL  CL?2D(DEBYE,SK,W,CHD,VP,NTOT,NBDW,NEL.NBDY,MBDY,XX,YY,MCOM, 
+ANEL,£NOT) 

C  READ(ll)  NXX.NYY, ( (WW(I,J) ,I-1,NXX) ,J-l,NYY) 

C  IP(MTHD.3Q.l)  THEN 

S  S-?30 

C  DO  130  I1-1,NX 

C  KX«KX-hl 

C  IF(NX.LE.NXX)  THEN 

C  WW(II,JJ)«WW{II,JJ)-VV(KK) 

C  W(KX)-VW(II,JJ) 

C  WW(II,JJ)-W(XK) 

C  SNDIP 

Cl 30  CONTINUE 
C  ENDIF 

C 

C  DO  CONTOURS 
NLB-12 
FCON(l)-VOUT 
FCON(2)— 1. 

PCON(3)-l. 

FCOH14^«VIHN 

IN(1)-1 

IN(2)-6 

IN(3)-7 

IN(4)-12 

CALL  VEXP(4,IN,FCON) 

IFPLT-0 

C 

C  DO  QUADRATURES 
IX-IAP(l) 
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o  n 


MMR»  ‘  MR-:  '  /C-^l 
MMA»‘  MA-:  :/ 

•MAT-: 

I?'M130  ..'TE.  :  '  THEN 
MAT=MA 
MMA-i 
MA=: 

ENDir 

MB-MAXf  MA.MAT', 

DC  30  D»1,MR 

R-R3/U 

DC  90  3  )  .lA?'-;) 

.W- 2 »  (  MX- 1 )  *  ( JJ- 1 ) -^2*  ( IX- 1 ) +1 

CALL  ACRDS(0.D0,XX(IX)  ,R.MN,XX,irY, W,MCON.AIJr/,AREA.ACHK) 
I?(AaS(ACHK-l. ) .-T.l.D-6)  GO  TO  110 
90  CCNTIMUE 

110  CCNTIMUE 

I? (MR  EQ.l)  UIMC«U*y/RB 
IE{MB.DE.l)  ALEIMC-U/RB 

DEPINITICNS  OF  ALPHA  QUADRATURE  PARAMETERS 
DR=-A3S(RB-R) 

CALF- ( R*R+RB*RB-RA*RA) /2 . /RB/R 
AL?A1--AC0S(MIM(1. ,CALF) ) 

ALFA2--.ALFA1 

DO  30  MM-1,MA 

IF(MiaO.SQ.2)  THEN 

ALFA-ALFAST+ ( MM- I ) *ALFINC 

I?(ALFA.LT.ALFA1.0R.ALPA.GT.ALFA2)  GO  TO  30 
SQ-SQRT(RC*RC-(RB*SIN(ALPA) )**2) 

RCAL-RB*CCS ( ALFA ) 

RIA-RCAL-SQ 
ROA-RCAL+SQ 
UIA-RB/RIA 
U0A-R3/R0A 
Ejroi? 

DO  40  I-l.ME 

EK- ( EKST+ ( I-l ) *EXINC) *SNOT 
DO  40  J-1,MP 

PHI-?HlST+( J-1 ) *PHINC 
XP-XX(IAP(1) ) 

YP-R 

MOWEL-NN 
NCWSD-2 

CALL  ACRDS(0. DO, XP,YP, MOWED, XX, YY.W.MCOM.AINV, AREA, ACHK) 
SPEED-3CRT(EX+AINV{3,5) ) 

XDOT»SPEED*COS ( PHI ) 

YDOT-S?EED»SIM( PHI ) 

MTRJ-NTRJ+1 

CALL  TRAJ  ( AINV , 0 . DO , XSTSP , XSGM , YSGM , W , XX . YY , IPRT , MTRJ , MEL 
♦  ,MC0M,:TEST, LOG, TAU,IFPLT,RB,AEND, SEND, BET. REND) 

DO  60  M-1,MAT 

:F(M130.ME.2)  THEN 

ALFA-AL?AST+(M-1 ) »ALFINC 
END  IF 

:F(MR. EQ.l. AND. MB. EQ.l)  THEN 
Il-I 


ELSE  ZE'ME.NE. 1 )  THEN 

r  1  *  r 

:?(MP  ..'IE.  1 ;  THEN 

ELSE  rE'MR.:iE.l)  THEN 
I2»L 
ELSE 

r2»MAX(M.MM) 

END  IF 

ELSE  IF(MP.NE.l)  THEN 
I2»J 

IFfMR.NE.l)  THEN 
I1«L 
ELSE 

END  I? 

ELSE 

11- L 

12- MAX(M,MM) 

MMR-MR 

ENDIF 

TEST( 11,12 )-'Q' 

r?(M130.NE.2.AND. (ALrA.LT.ALFA1.0R.AI.FA.GT.ALFA2) )  GO  TO 
+  60 

TEST (11.12) -ICHAR( ITEST ) 

IF( ITEST.NE. 1 )  GO  TO  60 

ISUM(I)-ISUM(I)+1 

JSUM( J)«JSUM( J)+l 

LSUM(L)-LSUM(L)+1 

MSUM(M)-MSUM(M)-<-l 

ASUM-ASUM+1 . /U/U/U 

ITOT  { I .  J ) -ITOT  ( I .  J ) -t-i 

NRHQ- ( REND-RI? ) / ( ROP-RIP ) »21 . + 1 . 

NRH0-M1M( MAX  ( 1  .NRHO ) ,  21) 

NSUM ( NRHO ) -NSUM ( NRHO ) + 1 

NEO-  ( E]END- 1000 .  -VOUT )  /  ( VINN-VCUT )  *21 .  +1 . 

NEO-MIN ( MAX ( 1 . NEO ), 21 ) 

NAO-(BET+8. )/16.*21.+l. 

NAO-MIN ( MAX ( 1 , NAO ) , 21 ) 

NEOT ( NEO ) -NEOT ( NEO ) +1 
NAOT ( NAO ) -NAOT ( NAO ) +1 
60  CONTINUE 

40  CONTINUE 

IF(L.SQ.MMR.AND.MM.SQ.MMA)  THEN 
IF(Il.GT.l.AND.I2.GT.l)  THEN 
DO  120  JJ-l.MIN(41.I2) 

WRITE(6,  '  (1X.41A3)  '  )  (TESTd.JJ)  ,1-1  ,MIN(  41 ,11 )  ) 

120  CONTINUE 

ENDIF 
ENDIF 

30  CONTINUE 

IF(Il.GT.l.AND.I2.GT.l)  THEN 
WRITE(6,33)  PHIl/RAD 

33  format;//'  phi  start-' ,F5.1, IX, 'DEGRESS' ) 

DO  30  :i-l.MIN(41,MP) 

WRITS(6, ' (IX, 4113) • )  (ITOT(I,II) ,I-1,MIN(41,ME) ) 

80  CONTINUE 

WRITE (6, 34)  PHI2/RAD 

34  format; IX, 'PHI  END-' ,F4.:, IX. 'DEGRESS' ) 
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6  ,  35  i  iX;,2K2 

35  rCRMAT;  IX,  'S.'fERGV  START- '  ,  J’ .  3 , 3SX  ,  '  SNEHGY  SND-',?7.5) 

end:? 

•^:te;  5 . :to  ) 

100  rCRMAT;//8X. 'E3EG‘  ,i:x, '?HI'  ,14X, 'R'  ,12X, ' AL? ’  ,13.X, 'SEND'  ,:iX, 
-  ■  3E.'rD  '  ,  :2X  ,  'RKC  '  ' 

DO  24  1 .MAX, MR.MB -ME.MP , 21 ) 

ITMP ( 1 , -  '  ' 

LTMP  f  2 1  - '  ' 

ITMP ' 3  ^ ' 

ITMP ; 4 1 > ’  • 

LTMP  '  S  ')  =  '  ' 

LTMP(6)-’  ■ 

LTMP ; 7 ) -  '  ' 

EX-EXST-^  '  )  *EXLNC 

?h:-?hist-^{  :-i )  *?htnc 

DPHI-PHT/RAD 

a-usT*(:-i)*u:Nc 

R-l  ./'J 


4001 

24 

6 

C 


23 


102 

9 

2001 


2002 


AL? A-ALFAST4-  ( I-l )  *AL?rNC 
SC-E3-*- FLOAT  ( I  - 1 )  *OE 
30--  3  .  -r  (  FLOAT  ( I )  -  .  S  )  *16  .  /21 . 

RO-RIP* I  FLOAT ( I ) - . S ) * (HOP-RIP ) /21 . 

:ss- r SUM ( I ) + JSUM ( I ) f LSUM ( I ) +MSUM  (I) 
rF(I.LE.  21 )  LSS-ISS+^^EOT(I)+MAOT(I)+^ISUM(I) 

IFCSS.LE.O)  GO  TO  24 
LFC.LE.ME)  WRITS(LTMP(1)  ,4001) 

IFC.LE.MP)  WRITS(LTMP(2)  ,4001) 

:F(I.LS.MR)  WRITS(LTMP(3) ,4001) 

IFC.LE.MB)  WRITE(LTMP(4)  ,4001) 

:?(I.LS.21)  WRITE(LTMP{S) ,4001) 

I?(:.LS.21)  WRITE(LTMP(6) ,4001) 
rF(I.LE.21)  WRITS(LTMP(7) ,4001) 

?CRMAT(I4,?10.S ) 

WRITE (6, 6)  LTMP 
CONTINUE 

FORMAT! 7 ( IX, A1 4) ) 

WRITE! 6, 101) 

IISM-0 

DC  23  1-1,11 

irsM-:isM+isuM!i) 

CONTINUE 

GF-  !ASUM)  •RB*RB*ST*ST»PHINC*UINC*SXINC*ALFINC 
THETA-THETA/ RAD 
WRITE (6, 102) 

F0RMAT!1X/1X/1X, •**»THETA**SUM  OVER  ALL  GEOMETRIC  FACTOR') 

WRITE ! 6 , 9 ) THETA , IISM , GF 

FCRMAT!1X,F10.S,I10,1?E15.5) 

IF!ME.GT.2)  THEN 

CALL  ?WKM!ME.ISUM,EXST.SXINC,SBAR.SIGE,ERAN,MAXX) 

SRAN2- FLOAT ( IISM ) *EKINC/FLOAT ! MAXX ) 

FORMAT! //IX, 'EBAR.SIGE.  DELTA  E  OVER  S  -  ’.4F10.S) 
WRITE! 6, 2001)  EBAR,SIGS.SRAN/MAX(EBAR,l.D-6) ,EHAN2 
SNDI? 

I?!MP.GT.2)  THEN 

CALL  FWHM !  MP ,  JSUM ,  PHIST ,  PHINC .  P3AR .  SIG? ,  PRAN , MAXX ) 
PRAV2- FLOAT !  IISM )  *PHINC/FLOAT  ( MAXX)  /RAD 
WRITE  (6,2002)  P3AR/RAD ,  SIG? /RAD ,  PSAN/RAD ,  PRAN2 
FORMAT ! //IX, '?3AR,3IG?,?FWHM  -  ',4F10.5) 

END  IF 


ISUM(I) ,2K 
JSUM(I) ,DPHI 
LSUM(I) ,R 
MSUM(I) , ALFA/RAD 
NEOT ! I ) , SO 
NAOT(I) ,30 
MSUM(I) ,30 
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;r;MR.3T.2)  then 

CAEE  ?WHM  (  MR  ,  13UM  ,  TST  ,  JINC ,  'JBAR  ,  SIGR  .  JRAN  .  MAXX ) 
•JRAN:-?L0AT  '  :iSM )  -UINC/FLOAT  (  MAXX 1 
•/raiTS,'  6  .  ;00  3  ')  1  ./UBAR.aTGR,  JRAN/JRAN2 
2003  FORMAT IX, 'RBAR,3IGR,RFWHM  -  ',4F10.5) 

FNDXF 

XF'MB.GT.:)  THEN 


CAL3  rWKM ( MB , MSUM , ALFAST , ALFXNC .  .\BAR , SIGA , ARAN , MAXX ) 

ARAN2-FLOAT  ( IISM 1  *AI.FTNC,'FLOAT  f  MAXX /RAD 
writs; S  ,2004)  AflAR/RAD,SIGA/RAD,ARAN/RAD,ARAN2 

2004  FCRMAT(//1X. ' ABAR,3IGA,AFWHM  ■  ',4F10.5) 

END  IF 

DE». 22/21. 

SB-.94*DE/2. 

CALG  FWHM( 21 .NEOT.EB.DE.EOBR.SIGE.SORAN.MAXX) 

WRITE;6,20QS)  ECBR.SIGS, SCRAN 

2005  FORMAT ( //IX. 'SCBR,SIGE,SFWHM  ■  '.3F10.S) 

DB-IS ./21. 

3B--3 . rDB/2. 

CALL  ?WHM(21,NACT.3B,DB.30BR,SIGa,B0RAN,MAXX) 

WRITS (6,2006)  BCBR.SIGB.BCRAN 

2006  FORMAT ( //IX, 'BOBR.SIGB.aFWHM  -  '.SFIO.S) 

DR»{ROP-aiP)/21. 

RHOB-RIP-*-DR/2. 

CALL  ?WHM(21,NSUM.RHOB,DR,aOBR,SlGR,aORAN,MAXX) 

WRITS(6,2007)  R0BR,SIGR,RCRAN 

2007  FORMAT ( //IX, 'ROBR,SlGR.aFWHM  ■  ',3P10.S) 

END 

SUBROUTINE  CNEC ( NX , NY , NEL , NBD Y , NBDW , MCCN , MBDY , ANEL , XX , YY , W ,  VP ) 
IMPLICIT  REAL'S  (A-H,0-2) 

DIMENSION  MCON(7,*) ,ANEL(») ,XX(«) ,YY(*) ,W(*) ,AINV(3,5) ,MBDY(2,*) , 
-VP ( * ) 

NEL-0 

NBDW-NX+2 

NTCT-NX'NY 

Nl»NTOT-NX 

IMlD-(NX-l)/2+l 

DO  10  NN-1,N1 

IN-MOD  ( NN- 1  ,NX )  -t-1 
IFdN.SQ.NX)  GO  TO  10 
NEL-NEL+2 
MCCN (1, NEL- 1)-NN 
MCCN  (  3 , NEL- 1 ) -NN+NX 
IF(IN.LT.IMID)  THEN 
MCCN ( 2 , NEL- I ) -NN+l+NX 
MCCN (1, NEL )-NN 
MCON(2,NEL)-NN+l 
MCCN( 3 ,NEL ) -NN+I+NX 
ELSE 

MCCN ( 2 , NEL- 1 ) -NN+1 
MCON(l,NEL)-NN+l 
MCCN ( 2 , NEL ) -NN+ 1 +NX 
MCCN ( 3, NEL) -NN+NX 
ENDI? 

10  CCNTINUE 
C 

DO  20  11-4,7 
DO  20  JJ-1,NEL 
MCCN(II,JJ)--1 
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o  o  o  o  o 


20 


CCNTINUE 


30  7  0  ::»L..'rEL 
30  ’0  JJ-L.j 

:?'MC3N(JJ*3,I3: .Or.O)  30  TO  70 
30  60  :<K-2,3TEL 

IF'XK.EQ.II )  GO  TO  50 
30  30  3L» I , 3 

:?  t  MOON  (  33  ,  :<X  )  .£Q.MCON(JJ.ri)  )  THEN 

POINT  33,;<X  -  JT.II 

30  40  MM»2,3 

go  to  40 

30  30  NN-1.3 

IJ(NN.EQ,33)  GO  TO  30 
IF(MCCN(NN,XK) . 2Q .MCON( MM , 11 ) )  THEN 

NN.XX  -  MM,  1 1  HENCE  XX  AND  II  ARE  ADJACENT  TRIANGLES  WITH 
SIDE  33-NN  OF  XX  --  SIDE  JJ-MM  OP  II 

11- MIN( JJ,MM) 

12- MAX( JJ,MM) 

X1-MIN(L3.NN) 

X2-MAX(L3.NN) 

IS-3-I1 

ir(IS.SQ.I2)  IS-3 
XS-3-X1 

I?(XS.EQ.X2)  XS-3 
MCON(IS+3,II)-KX 
MCQN(XS+3,XX)-II 
ENDI? 

30  CONTINUE 

40  CC''^INUE 

SNDIF 

SO  CONTINUE 

60  CONTINUE 

70  CONTINUE 
C 

DO  30  II»1,NE3 
DO  80  JJ-4,6 

:P(MC0N( JJ. II) .GT.O)  GO  TO  30 

11*1 

12"  3 

IF(JJ.EQ.4)  11-2 
IF{JJ.S0.6)  12-2 
J1-M0D(MC0N(I1,II)-1,NX)+1 
J2-MOD(MCON(I2,II)-l,NX)-t.l 

I?(MC0N(I1.II) .3E.NX.AND.MCON(I2,II) .3E.NX)  THEN 
C  LOWER  BOUNDARY 

MCON(JJ,II)-0 

ELSE  I?(MCCN(I1,II)  .GT.NTOT-NX.AND.MCON(I2,II) .GT.NTOT-NX)  T 
+  HEN 

C  UPPER  BOUNDARY 

MCON( JJ,II)--1 

ELSE  I?{Jl.S<3.1.AND.J2.SQ.l)  THEN 
C  LEFT  END 

MCCN(JJ,II)--2 

ELSE  IF{J1 .EQ.NX.AND.J2.EQ.NX)  THEN 
C  BIGHT  END 

MCGNf  JJ,:i>-? 

ENDI? 
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30  cc>rrrNU£ 

C  CONNECT  ~rXED  POTENTIALS  EHCM  MCON 
00  ?0  ;i«;,NBDY-l 
:?T=MBDTff  i,ii; 

OC  100  TJ-:,NEL 
00  100  ;<K=«1,3 

IF'MCCNtXK.JJ) -EC- -?T',  THEN 
:?»M0D(KK.3)-1 
T?»MOD(IP,3)-rl 
00  120  LL-II-1,NBDT 

:F(MBDY(2.1L) .SQ.MBDT(2,II) )  THEN 

IFfMCOMdP.OJ)  .SQ.MBDY(1,LL)  )  MCON(  JP+3  .  JJ) -MBDY  (  2  ,LL ) 
IF(MCON(OP,JJ)  .SQ.MBDYd.LL)  )  MCCN(  IP-t-3  ,  JJ) -MBDY  (  2  .LL ) 
SNDIF 


120 

CONTINUE 

GO  TO  100 

ENDIF 

110 

CONTINUE 

100 

CONTINUE 

90 

CONTINUE 

:  REMOVE  -99.  MARKED  FIXED 

”■0 


00  140  II-1,NBDY 

IF(VP(II) .NE.-99. )  THEN 
JJ* JJ+l 

MBDY ( 1 , Jw ) -MBDY (1,11) 

MBDY ( 2 . JJ ) -MBDY ( 2 . II ) 

VP(JJ)-VP(II) 

ENDIF 

140  CONTINUE 
NBDY-JJ 
C 

C  CALCULATE  AREAS  OP  TRIANGLES 
DO  130  II-1,NEL 

CALL  ACRDS(O.DO,O.DO,O.DO,I1,XX,YY,W,MCON,AINV,ANEL(1I) ,ASUM) 
130  CONTINUE 
C 

RETURN 

END 

SUBROUTINE  ACRDS ( VPRP ,XPT , YPT ,NIE , XX , YY , W ,MCON.BINV,AREA.ASaM) 
TMor  pp&r*fl  ^ 

DIMENSION  XX(*) ,YY(*) ,AINV(3,5) .BINV(3,5) ,MCON(7,*) . 

•rXEC  3  )  ,YE(3)  ,Wf*) 

DOUBLE  PRECISION  A2,XE,YE,AINV 
C 

ZRO-l.S-9 
n  10  1-1,3 

XE(I)-XX(MC0N(I,NIE) ) 

YS(I)-YY(MCON(I,NIE) ) 

10  CONTINUE 
C 

AI NV  ( 1 , 1 )  -  XE  (  2 )  *  YE  (  3  )  -  XE  (  3  )  *  YE  ( 2 ) 

AI NV  (  2 , 1 )  -  XE  (  3  )  *  YE  ( 1 )  -  XE  ( 1 )  *  YE  (  3  ) 

AINV (3,1)»XE(1)*YE(2)-XE(2) »YE ( 1 ) 

A2- AINV (1.1) +AINV (2.1) +AINV (3.1) 

AREA-A2/2. 

AINV (1,1) -AINV ( 1 . 1 ) /A2 
AINV (2.1) -AINV ( 2 , 1 ) /A2 
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A:>rv  »AINV  (3,1:  /A2 

AINV  f 1 , 1 ! -  I YS  f  2 1 -  'S ( 3 1  ) /A2 
A2>rv  {  2 , 2  1  »  '  YS  t  3  )  -  Y£  (  1 )  )  /A2 
Al.'rV  '  3  .  2  !  -  '  YE  f  1  -YE ;  2  )  1  /A2 
A2.W  f  1  ,  3  1  »  '  XE  (  3  )  -:<E  (  2  ',  )  /A2 
A2NV  r  2 , 3  )  =  '  XE  ^  1 ;  -  :<E  (  3  1  1  / A2 
AI.W  (  3  ,  3  1  =  I  ;<E 2  :  -XE  ( 1 )  )  /A2 

AIXV'l,3i>iO. 

A2NV( 2 , 5  i-O . 

A2.'rV(  3 , 5  (-0  . 

20  20  :=i,3 

AINV  tl.i)  -AINV  (2,1)  +AINV  (1,2)  «X?T-‘-AINV  (2,3)  ♦YPT 
20  20  J«2,4 

22-MC0N(2,.'r:E) 

.A2.W  ;  J-1 ,  S  )  -A23fV{  J-l  ,5  )  +A2NV(  2  ,  J  )  ♦  ('/V(  II )  -r  ( VPSP/MAX(  .  001  ,ABS  ( 
YY(22) ; ) )**2) 

20  CCNI'INOE 

AS'JM-0  . 

30  30  2»1,3 

2?(A3S(A2>r7(I,4  )  )  .31.250)  AINV  ( 1 , 4  )  *0  . 

2?(ABS(AINV(I.5) ) .3T.2S0)  AINV(I,5)-0. 

AS1IM-ASUM+ABS(AINV(I,4)  ) 

30  CONTINUE 

DO  40  22-1,5 
DO  40  JJ-1,3 

3INV{JJ,II)-A2NV(JJ,II) 

40  CONTINUE 
C 

RETURN 

END 

SUBROUTINE  CL?20{  DEBYE ,  SK ,  W ,  CHD ,  VP  ,  NTOT ,  NBDW ,  NEL ,  NBDY , MBDY ,  XX ,  YY , 
tMCON.ANEL.ENOT) 

IMPLICIT  REAL*a  (A-H,0-Z) 

DIMENSION  SK(NBDW,NTOT)  ,VV(NT(0T)  ,CHD(NEI.)  ,VP(NBDY)  ,XX(NT0T) , 
^YY(NTOT)  ,ANEL(NEL)  ..MBDY (2, NBDY)  ,MCON(7,NEL) 

T’/OLTS-AaS(ENOT) 

IEAXIS-3 

DO  10  II- 1,. NBDY 

V? ( II ) -VP ( II ) /TVOLTS 
10  CONTINUE 
C 

C  SET  RIGHT  HAND  SIDE 

CAL3  3ETQ2(NT0T,NEL,DEBYE,VV,CHD,ANEL,MC0N) 

C 

C  SST  SK 

*CAL3  SETX0( NTOT, NBDW, NEL, IEAXIS,.MCCN,SX, XX, YY,ANEL) 

C 

C  IMPOSE  ESSENTIAL  B.C. 

CALL  ESSaC  ( SX  ,  W ,  VP  , NBDY , .MBDY , .NTOT ,  .NBDW ) 

C 

C  SCL’/E 

CALL  SOLV(SX,V'/,. NTOT,. NBDW) 

C 

C  CONVERT  FROM  XT  TO  VOLTS 
DO  20  I- 1,. NTOT 

IF'ABSC/VC)  )  .LT.l.E-9)  VV(I,-0. 

■/V  ( I )  -7V  { I )  *TVCLTS 
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o  o  o  o  o  no 


20 


"TVOLTS 


r? ' : .  IE .  )  vp « I )  «vp  >  2 ; 

CONTINUE 
EETURN 
END 

SUBROUTINE  TRAJ  !  3INV  ,  7PRP  ,  XSTS?  ,  XSGN .  ISGN ,  'ri .  XX ,  Ti  ,  I?RT .  NTRJ  , 

-NEi. , MOON , lEND , IOC . TAU , lEPLT , RB , AEND , SGTS , 3ETS  ,RPS ) 

IMPLICIT  REAL*9  (A-H,0-2) 

DIMENSION  MC0N(7,*)  ,XX(*)  ,TY(*) ,TAUf*) ,AINV(3,3) , 

-?LT(2,51) ,3INV(3.S) 

COMMON  /ORB/  XDOT  ,  TDOT . XP  ,  TP  , NTSAV  ,  TIMIN , NOWEL , NOWSO 
COMMON  /INST/  RIP,ROP.aiA.RCA.ENOT,3G2 
SAVE 

DATA  ICALO  /O/ 

I?(IC\LL.EQ.O)  THEN 
ICALL-l 
?I=«ACOS(-l.  ) 

RTCD-iaO./PI 
C  CALL  SSWTCH  (  2  ,  IS'W2  ) 

ISW2»2 

RKON- ( 1 . /RIP-1 . /ROP ) 

EXON- SNOT* ( ROP+RIP ) *RXON 
VINN-EKON/ ( 1 . +RIP/RCP ) 

VOUT-VINN-SKON 

ENDIF 

C 

?LT(l,51)-0. 

?LT(2,Sl)-0. 

NPT-0 

:<STSP-o 

XSMX-IOOO 

TSGN-1. 

XSGN-1. 

ERMX-0. 

RP-SQRT ( XP*XP+TP*TP ) 

ALP-ATAN2 ( TP , XP ) *RTOD 
I?(AL?.GT.1S0. )  ALP-ALP-360. 

BET- 9  0 . +ATAN2 ( TDOT , XDOT ) *RTOD- ALP 
IP(BET.GT.90. )  BET-BET-360. 

IF(BET.LT.-90. )  BET-BET+360. 

EGT»XDOT*XDOT+Tr)OT*TDOT 
I?(IPRT..NE.O)  THEN 
WRITE (6, 3002) 

WRITE(6,300l)  KSTSP,.NOWEL.NCWSD,XP,TP,XDOT,TTX3T.BINV(3,S)  ,IEND, 
0,0.0.  .RP, ALP, EGY, BET 

3001  FORMAT  ( IX  ,  14 , 14 , 13 , 2F8 . 4 , 2F3 . 2  ,  F8 . 4 . 14 , 214  .F8 . 4  ,F10 . 4 , 2F8 . 2  ,F8 . 4 

3002  FORMAT(//3X, 'KS  NIE  NSD  X  T’,6X.'  ,4X,'YDOT', 

+  5X, 'PHI' .2X, 'IE  IND  LOG  T' .9X, 'R' ■ ,5X. 'EGT* ,4X. 'BET* 

+  ) 

ENDIF 

CHECK  FOR  REFLECTIVE  INJECTION 
INDSV-1 
LOC-0 

.NOT  NOW  USING  T-0  SYMMETRY  PLANE 

IF(TDOT.GE.O.  .OR.YY..NE.O.  )  GO  TO  60 
TSGN— 1. 

TDOT— TDOT 

60  CONTINUE 


0(10  o  o 


.'fCT  ::cw  jsi.vG  :<=•  :  symmeth^;  plane 

:E'XC0T.3E.3. .CR.XP.NE.O. )  GOTO  10 

xsg::=.-i  . 

XLCT  =  -Xi;CT 

10  GCMTl.'fUE 

.'rE:<EL=MCCN  (  NCWSD*  :  ,  N'CWEL 
:r '  KSTE? .  irE .  0 )  then 
3C  30  :s»l,3 

:E  '  MCON  (  rs-r  3  ,  NEXEL  )  .  eq  .  ncwel  >  ncwsd-  is 
30  continue 

NCWEL=NEXEL 

Eiroi? 

20  CONTINUE 
IERH=0 

CALL  CR3 :T2  AIN%' ,  I?RT  ,  VPRP ,  XSTSP  ,  MCON ,  NEL  ,  W  ,  XX ,  YY  . NTRJ ,  lEND  , 
•*LCC,ERMX,3I>r7) 

;<STE?  =  ;<STEP*1 
I?(I??LT.NE.O)  THEN 
NPT  •  NPT  1 
?LT( 1 ,NFT) -XP 
?LT( 2,NPT)»YP 
I?;i?PLT.EQ.2)  THEN 

PLT ( 1 , NPT ) - ATAN2 ( Y? . XP ) *RTCD 
PLT ( 2 , NPT ) -SCRT ( X?*X?+YP*YP ) 

I?'?LT( 1 ,NFT) .GT.90. )  ?LT( 1 ,NPT ) «?LT ( 1 ,NPT ) -360 . 

ENCI? 

ie:npt.sq.50)  then 
WRITEO)  2, 51,  PLT 
.'rpT»l 

?LT'1,1)-?LT{1,50) 

.';2,1)«?LT'2,30) 

”...  •  ,j  I E 
END:- 

INi:iC*MC0N(NCWS0*3  .NCWEL) 

C?RT*0 

I?; IPRT.EC. 1 .OR. (IPRT.GT.I.AND.INDIC.LE.O) )  JPRT-1 
I?'AaS(XP) .1T.1.3-6.0R.ABS(YP) . LT . I . E-6 .OR. INDIC.LE. 0 )  THEN 
RP>SCRT!X?*X?*YP*YP) 

AL?^ATAN2 ( Y? ,XP ) *HTOD 
lEfALP.GT.lSO. )  ALP-ALP-360. 

3ET  -  90  .  -^ATAN2  ( YDOT ,  XDOT)  *RTOD-AL? 

I?(3ST.GT.90. )  3ET-3ET-360. 

IFf3ET.LT.-90. )  3ET-3ET+360. 

EG Y - XDOT »  XDOT + YDOT* YDOT 
I?'A3S(AL?-AEND) .LT.l.E-5)  THEN 
RPS-RP 
EGYE-EGY 
3ETE-3ET 
ElTDIE 

I?'I?RT.EQ.2)  JPRT-1 
ENDir 

I?{ JPRT.EQ. 1 )  THEN 

WRITS(  6  ,  3001 )  XSTSP, NCWEL,.NCWSD,X?*XSGN,Y?»YSGN,XDOT*XSGN, YDOT 
-  *YSGN,3I.NV(  3 ,3)  ,IEND,INDIC,LCC,TIMIN,SP,ALP  "IGY.BET 
Eire  IE 

TAU  '  .VOWEL  )  -TAU  f  NOWEL )  +TIMIN 
lE'LCC.EO.O)  I.VDSV=1 
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liND .iQ . -1 )  GO  TO  90 
17 ( TEND. IT. 0)  GO  TO  110 
:7(KSTSP .EQ.KSMX)  GO  TO  SO 
17' INDIC.EQ. -1 )  GO  TO  100 
17':XD1C.GT.0)  GO  TO  10 

lE'lCC.N’E.O.ANC.  ;  I'rOIC  10  •  -  •  TR .  l.VDIC .  EQ  .  -  2  )  )  THEN 
17':.NTirC.EQ.rNDSV)  GO  TO  90 
INDSV.INDIC 
.'JCWSD-NTSAV 
GO  TO  20 
ENDI7 

C  END -OF -TRAJECTORY 
lEND-INDIC 
GO  TO  110 
C 

100  CONTINUE 

IF(A3S(XP) .lE.l.E-S)  THEN 
C  Y-AXIS  REFLECTION 
XSGN--XSGN 
XDCT— XDOT 

ELSE  I?{AaS(Y?) .LS.l.E-6)  THEN 
C  X-.\XIS  REFLECTION 
YSGN--YSGN 
YDCT=i-YDOT 
ELSE 

GO  TO  40 
ENDI7 

NCWSD-NTSAV 
GO  TO  20 
C 

so  CONTINUE 

C  TCO  MANY  STEPS.  TREAT  AS  ABSORPTION. 

C  WRITS (7,*)  •  ORBIT  STEPS  SCSTEP.XSMX  -  '  ,XSTEP,KSMX 

lEND— 3 
GO  TO  110 
C 

40  CONTINUE 

C  WRITS(7,*)  '  ORBIT  ERROR.  NTRJ  -  ' ,NTRJ 
IEND--3 
GO  TO  110 
C 

C  INJECTION  INTO  BODY.  TREAT  AS  ABSORPTION 
90  CONTINUE 
IEND--2 

IF(INDSV.LE.O)  lEND-INDSV 
C 

110  CONTINUE 

IFf IEND.EQ.-9)  THEN 
JALP-ALP 
DO  120  11-1,3 
XSTSP-KSTEP-vl 

CALL  ELLPS ( SGY , RP , BET . ALP , DAL? , SGYE . RPS , 3ETE . ALPS , LTEST ) 
INDrC--LTEST 

C  WRITE(6,3001)  XSTSP,. VOWEL, NCWSD,X?*XSGN,YP*YSGN,XDOT*XSGN,YDOT 

C  *YSGN,3INV{3,5)  ,IEND.INDIC,LOC,TIMIN,RPS,AL?E,EGYE,BETE 

ALP-ALPS 
DALP-90. 

RP-RPS 


86 


o  n  o  o  n 


130 


130 


C 

c 

c 


c 


SGY « cGYS 
3ET-3ETH: 

:r'GTESr..VE.:  ) 

GGiVTINL’S 
GG>rrrNUE 
E.'rarr 

:sw:.rc. : .  vtrite.'  5 .  jOod  .<srsp,:'iowEL,.'JcwsD,xp»xsGN, yp*ysgn, 

-XGCT^XSGIJ.  YECT'*YSGN,3XNV(  3,5)  ,  ISND  ,  INDIC  ,  IOC  ,  TIMIN 

:e  NPT.GT.i;  wRrrE;3)  x.npt-^:,  i  (Plt(:,J)  ,:=i,2)  ,j=«i,NP'r)  ,o.  ,o. 

MTEST=: . 

XE'AL?-:  .Z-3  .GT.AENDl  "THEN 

:E'AaSfAL?-9Q. ) THEN 
MTEST»4 

EOoE  :?'AL?.GT.-90. )  THEN 
MTEST=  2 
ELSE 

yiTSST=-6 
END  IE 

IE(H?.LT. {RI?*RCP)/2. )  MTEST=MTEST+1 
ENDIE 
lEiro-MTEST 
RETURN 
Eiro 

SUBRCCTI NE  CRB IT 2 ( AINV , IPRT , VPRP , XSTEP , MCON , NEL , W , XX , EY , NTRJ , 
*;<ESR  ,  ICC  ,  ERMX  ,  3INV ) 

IMPLICIT  REAL*3  (A-H.C-2) 

DETERMINATION  CF  PHIX,  PHIY ,  AT  X,  Y  FOR  NEXEL 

CC^^MCN  /CRB/  XECT,  YDOT,X?,Y?.NTSAV,TIMIN,NOWEL,NOWSD 
DOUBLE  PRECISION  ROOT, RADIO, TM,UD,UDD, ADD, D,D2RO,T, TIMM, XD, YD, 
-XNEW , YNEW , DY? , DX? 

DIMENSION  MCON(7,*) .'Ar(*) ,XX(*1 .YY(*1 
DIMENSION  T(6)  ,D(3)  ,’JD(  3  )  ,UDD(  3  ) 

DIMENSION  AINV( 3,5) ,3INV(3,S) 

SAVE 

DATA  IER,ICALL  /2*0/ 
lEf ICALL.EQ.O)  THEN 
LCC-0 
ICAI,L-«  1 
ENDIE 


IE»NOWEL 
TIMIN-0 . 

TIMN-0 . 

CALL  ACRDS ( V?RP , X? , Y? , NCWEL , XX , YY , W , MCON , AINV , AREA , ASUM ) 

XD-XDCT 

YD -Y DOT 

DX?«X? 

DY?«YP 

XERR-0 

IALT-0 

NTSAV»NCWSD 


STEP  ACROSS  2-D  TRIANGLE  ASSUMING  CONSTANT  POTENTIAL  WITHIN  TRIANGLE 
INCLUDING  ACCELERATION  TERMS 


TCCM=3 . 3333E+33 
2RC»I .E-3 


n  o 


DZRC»: . 


;<DD=  .  5*A:.Wf  1,5  ', 

•fDD=.3*AI.'rv»;,5) 

c 

70  CCNTINUS 

00  10  :=»i,3 

0(1: =AINV( I , 4 ) 

■JD(  r  »AIMVr  I  ,  0)  *XDOT-^AINV(  1 , 2  )  *''OCT 
UDD  (  I  1  »  AINV  (1,0)  *XDD+AINV  (1,2)  *''DD 

10  coim.vi:s 

c 

ISIO-O 
DO  50  :<»1.2 

:?(0(X) .iT.oac)  D(X)-o. 

:?!AaS(UD(:<)  )  .LT.2R0)  UD(K)-0. 

I?(UDD(X)*UDD(X)  .DT.OZRO)  •JDD(K)-0. 

I?(D(X)  .EC.O,  .AND.K.JIE.NCWSD)  THEN 

:f(Ncwsd.ne.O)  isid-x 
I?(NCWSD.2Q.O)  NCWSD-K 
ENDI? 

50  CONTINUE 

I?(ISID.GT.O.AND.NCWSD.GT.O)  THEN 

I?!UD(NCWSD) .DT.O. .OR. (UD(NCWSD) .2Q.0. . AND . UDD ( NOWSD ) .LT.O. ) .OR 
+  UD(ISID) .LT.O. .OR. (UD(ISID) .EQ.O. .AND.UDD(ISID) .DT.O. ) )  THEN 
.NEXSD-MCON ( ISID+  3  . NCWEL ) 

IFdCC.GS.lO )  WRITE(7,«)  '  LOG,. NCWEL,  NCWSD,ISID,NEXEL,UD  -  ' 
■(■  ICC ,  .NCWEL ,  ISID , .NEXEL ,  UD ( NCWSD )  .  UD  ( ISID ) 

I?(LCC.LE.16)  THEN 
LCC»LCC-t-l 
NTSAV-’SID 
NCWSD- ISID 
ELSE 

XESE--1 

LCC-0 

SNDI? 

RETURN 

ENDI? 

ENDI? 

C 

LCC«0 

DO  20  X-1.3 

I?(0(X)  .NE.O.  .OR.UD(X)  .GT.O.  .OR.ISID..NE.O)  GOTO  20 

I?(XSTSP.EQ.O)  RETURN 

UD(X)-l.E-4 

I?(AINV(X,2) .EQ.O. )  GO  TO  30 
XD-(UD(K)-AINV(X,3)*YDOT)/AINV(X,2) 

XDOT-XD 
GO  TO  70 
C 

30  CONTI.NUE 

I?(AINV(X,3) .NE.O. )  YD-(UD(X)-AINV(X,2)*XDOT)/AINV(K,  3) 
I?(AINV(X,3) .NE.O. )  YDOT-YD 
20  CONTINUE 

I?(IER.EQ.l)  WRITE(6.300I)  D,UD 
3001  FORMAT ( SX, 'D.UD  -  ' , 2 ( 3F10 . 6 . 2X ) ) 

C 

NTSAV--9 

TIMN-TCCM 
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n  o  o  no 


3C  503  :<=«:,3 

x:»x2-: 

T'Kl  )-rCOM 
T  '  K2  -TCCM 
RCCT-0 . 

:?"JDD^:<)  .SQ.O.  .AND.2’DfK;  .ZQ.3.  )  30  TO  500 
;-''JDD(:<l  .:fE.3.  .CR.OD(K)  .ZQ.O.  )  GOTO  100 
TM-'OCX', /UD(K', 

GO  TO  200 

100  ZCNTIlfUZ 

rM--UD(  :<)  /UDDCO 
AOD--2  .  •*D(:<)/UDD(:<) 

RADIC-TM«TM-t-ADD 

;Z f RADIO. Lr.-02RC)  GO  TO  500 

RADIO.  GT.D2RO)  RCCT=»DSORT(  RADIO) 

200  00NTINU2 

IZ;:<.2Q.l)  WRITE(7,*)  •  TM, ADD, RADIO, ROOT, T  ■  ’, TM , ADD , RADIO, 

■r  RCCT,T(X1)  ,T(X2) 

T(X1 ) -TM-RCOT 
Tf X2 ) =TM+RCOT 
00  CONTIlfUE 

RIND  SHORTEST  TIME  AfTO  SAVE  OCORDS . 

DO  400  D»l,6 

.'rr«(D-i)/2+i 

I?;T(L) .LT.O. .OR.T(L) .GE.TIMN)  GO  TO  400 

I?(  (NT.EQ..VCWSD.QR.NT.EQ.ISID)  .AilD.T(D)  .LE.ZRO)  GO  TO  400 
0  WRITS(7,»)  'NT,NOWSD,ISID,TIMN,T(L)  -  * ,OT,NOWSD, ISID,TIMN,T(I,) 

.vtsav-nt 

timk»ti;l) 

400  CONTINUE 

TIMIN»TIMN 

lE'XERR.SQ.l)  WRITE(7,*)  '  T  »',T 
XNEW--1 . 

TNEW--1 . 

3SUM»-1 . 

XDCT--1 . 

YDCT=«-1 . 

C 

IF(NTSAV.DT.O)  GO  TO  700 

0 

C  advance  through  minimum  TIME  TIMN 
IERR*0 

XNEW-DX?+XD*TIMN-*- .  S  »XDD*TIMN*TIMN 
INEW-DYP-^'jrD*TIMN+ .  S*YDD*TIMN*TIMN 
XT-XNEW 
YT-INEW 

0 

0  CHECX  PRECISION  OP  FINAL  DESTINATION 
C 

CALL  AC3DS<0 .DO ,XT,IT ,NCWEL,XX , YY , W ,MCON , 3INV , AREA, 3SUM) 
ZRMX»MAX(ABS(3SUM)-1. ,ERMX) 

I?;3SUM.GT.l.-«-25.*2RO)  IERR-1 

0 

0  FINISH  UP 
0 
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:<DCT  =  :0 -t- XD  D  *  T I MN 

xdot=yc+';2d*timn 

IF' rHRR..'fE.  J  )  GO  rc  '00 

C  SET  J?  SCOKXEEPIMG  ?CR  ADVANCE  TC  JIE.XT  TRIANGLE 
NCWSGaNTSAV 
;<?  =  XNEW 
T?=INEW 

IF'ABS,  XFSLT.2RO)  X?.0  . 

IF'ABS.Y?) .LT.2RO)  T?-0. 

IF  { A3S  X2CT  )  .  LT  .  IRC  1  XDCT-0. 

IFfAaS(ETOT) .LT.2RO)  E^CT-O. 

C 

RETURN 

700  CONTINUE 
lER*'. 

WRITS( 6 ,2001)  IERR,I3.NCWS0,NTSAV,ISI0,1CC,NTRU 

2001  FORMAT  ( //5X, 'PROBLEM  IN  CR3IT2  'i^ITH  lERR,  IS,NCWSD,NTSAV,ISID,LCC  ■ 
+  '  ,513,  'NTRJ  «  ■ ,17) 

WRITS ( 6,2002)  XP , YP , XD , ED , XDD , EDD , XNEW , ENEW , XDCT , YDOT , D , UD , ODD , 
*(AINV(I,4) ,I-1,3),ASUM-1..{BINV(I,4) ,1-1,3) ,3SUM-1. 

2002  FORMAT (/SX, 'INITIAL  X , E , XDOT , EDOT , XDD , EDD  -  ■,S?12.6/ 

+5X,’  FINAL  X,E, XDOT, EDOT  -  ',4F12.S/ 


*5X,  'D  '/ECTOR 

-  ',1P3E12.4/ 

-5X,'UD  VECTOR 

-  ',3E12.4/ 

-5X,'UDD  VECTOR 

-  ',3E12.4/ 

-SX,'A  VECTOR 

-  ',3E12.4,' 

ASUM-1.  - 

' ,£12.4/ 

-r5X,'3  '/ECTOR 

-  ',3E12.4,' 

BSUM-1.  - 

' ,E12.4) 

IF ( NTSAV . GT . 0 ) 

'WRITS  (6, 200  3  ) 

T,TIMN 

IF(NTSAV.LE.O) 

WRITE (6, 200 3) 

T,-l. 

2003  FORMAT (SX,'T  VECTOR  -  ', 1P6E12 . 4/5X ,' SELECTED  TIMN  -  ’,B14.6) 
KERR- -3 

C  WRITS  re,*)  '  0R3IT2  ERROR,  TRAJECTORE  NO  •  ' ,NTRJ 
C  CALL  DISCON(M) 

RETURN 

END 

SUBROUTINE  GECMB  ( VOUT ,  VINN ,  RIP ,  ROP ,  RIA  ,SOA ,  ENOT , NX .  NY , NBDE ,  XX ,  EE , 
+ W , CHD , VP , MBDE , PP , RR , lA? , AEND , IPRT , I3K , NBK ) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  XX(*) ,EE(«)  //V(*) ,RR(*) ,IN(65) ,XA( 65 ) .XT ( 65 ) ,SA( 65 ) , 
t-ET(65)  ,CHD(*)  ,VP{*)  ,mBDE(2,*)  ,  ?P  (  •  )  ,  IAP  (  4 , 4  )  ,  PA(  65  )  ,IBK(*) 
CHARACTER* 80  IHED 
C 

RTOD-190./ACOS(-L. ) 

C 

C  READ  INSTRUMENT  HEADER  CARD 

READ(4, ' (A80) ' ,SND-L00)  IHED 
WRITE(6, ' (1H1//SX,A) • )  IHED 
C 

C  READ  LOGIC  CARD 

READ(4 ,200L,END-100)  NX,NY,NL,NV,NIA,NOA,NCA,I?RT 
2001  F0RMAT(4I5,2X,3I1,I5) 

C 

C  READ  INPUT  APERTURE  CARD 
READ(4,*, END-100)  ZAP 
C 

C  READ  R  CARD  INPUT  AND  EXPAND 

3EAD( 4 , * ,END-100 )  NRP , (IN(I) ,RA(I) ,I-1,NRP) 

CALL  VEXP(NRP,IN,RA) 

C 


90 


RIP'-RAl  ZAP!  2 .2)) 

RCP’fUf  lAPf  4 . 2  )  ) 
piA=*^(  :aP!  : , : ) ) 

?.CA»SAl  :AP  '4,11) 

•^ITS  ;  6  .  '  ■  1 2X  ,  '  '  DIAMETER  OF  ENTRANCE  PlPERTURE  «  '  '  ,  F7 . 3  ,  '  '  CM '  '  ) 
-RCA-RIA 


•WRITE;  6 

,  '  '  5X,  '  ' INNER 

.\ND  CUTER 

RADIUS 

OF 

INSTRUMENT  » 

’  ' , 2F7 . 3  , 

'  OM'  ' 

■write;  5 

)  '  ;  RIP.ROP 
,  '  (3X. ■  'INNER 

.AND  CUTER 

RADIUS 

•AT 

130  DEG  -  ' ’ , 

,  2F7 . 3 

- ! ' ,  RAi 

writs;  s 

lAP ' 3 , 3 ) ) ,RA( : 
,  '  ;  7X,  •  ’  I.NNER 

lAP (4,3)) 
AND  OUTER 

RADIUS 

AT 

EXIT  APP  -  ' 

'  ,2F7 . 3 

RA(rAP( 3,4)) .RA(IA?(4,4) ) 


C 

RKCN= ( 1 . /RIP-1 . /RCP ) 

SXCNaSNCT* (ROP*ai?) *RXCN 
VI.VN-EXCN/  ( 1 .  -t-RIP/ROP  ) 

VCUT-VINN-SXON 

C 

C  READ  ALPHA  INPUT  CARD  AND  EXPAND 

READ (4,*, END- 100)  NAP, (IN(I) ,PA(I) ,1-1, NAP) 

CALL  VEX? ' NAP , IN , PA ) 

AEND-?A(NX) 

WRITE ( 6, ' (19X, * ’ANGULAR  LENGTH  OF  TUBE  -  •■,F7.1,''  DEG’’)')  90. 
*AEND 
C 

C  READ  BASE  LE’/EL  CARD  AND  EXPAND 

READ(4,*,END-100)  .NX?,  (IN(I)  ,XA(I)  .I-1,NXP) 

CALL  VEX?(.NXP,IN,XA) 

C 

C  DO  ALPHA  EXPANSION  OF  XA 
RTOD-130./ACOS{-i. ) 

DO  10  II-1,.NX 

IF'PAdI)  .SQ.90.  .OR.PAdI)  .EQ.-90.  )  THEN 
XXdD-XAdl) 

TTdl)-RA(l) 

IF(PAdl)  .SQ.-90.  )  YY(II)--RA(1) 

RR { II ) -3CRT { XX (11 ) **2+rY  ( II ) **2 ) 

?PdI)-ATAN2(YY(II)  ,XX(II)  )  *RTCD 
ELSE 

XX (I I ) -RA (1 ) *COS ( PA ( II ) /RTOD ) 

YY ( II ) -RA ( 1 ) *SIN( PA ( II ) /RTOD ) 

RRdl)-RA(l) 

PPdD-PAdI) 

END  IF 

I?(PPdI)  .GT.ISO.  )  PP(II)-?P(II)-360. 

10  CONTI.NUE 
C 

C  READ  TOP  LEVEL  CARD  AND  EXPAND 

READ( 4,*, END-100)  NX? . dN( I ) ,XA( I ) , I-l ,NX? ) 

CALL  '/EXP(NXP,IN,XA) 

C 

C  DO  ALPHA  EXPANSION  OF  XA 
DO  20  II-1,NX 

I? (PA (ID  .SQ.90.  .OR.PAdI)  .EQ. -90.  )  THEN 
XTdl)-XA(II) 

YTdl)-RA(.NY) 

IF(PA(II) .EG.-90. )  YT(II)--RA(NY) 

ELSE 

XT  (II )  -RA  ( .NY )  »COS  (  ?A(  II )  /RTOD ) 
YTdI)-RA(NY)*SIN(?AdI)/RTOD) 
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0 


£Ni::~ 

CC.TTI.'fUE 
C 

C  F.X?.!VNt)  3RI2  ?RCM  3CTTCM  TO  TCP  LEVELS 
TC  30  ::»2,NY 

RTrC= I RAi :i) -RA( L ) ) / ( RA ( NY ) -3A ( 1 ) ) 

LC  40  TJ=2,NX 
NN  =  NN-<-- 

XX ;  NN  ^  »  XX  (  J J )  +RTrO*  '  XT  (  J  J  )  - :<X  I J J 1  ) 

TY ; NN 1 »  Ti  ( J J ) +RTro* ( YT ( J J ) - YY ( J J )  ) 

?? ( NN ) = ATAN2 ( YY ( NN ) i XX ( NN ) ) *RTCD 
Tr fABS(XX(NN) ) .LT.l.S-10)  XX(NN)-0. 

IJf ABS( YY(NN) ) .LT.l.S-10)  YY(NN)-0. 

IS ( ASS ( ?9 ( NN ) ) . LT . 1 . 3- 1 0 )  PP ( NN ) *0 . 

IS(?P'.NN)  .GT.150.  )  PP(NN)-?P(NN)-360  . 
RR(.NN)a3CRT{XX(NN)**2+YY(NN)**2) 

40  CONTT.NUE 

30  CCNTI.NUE 
C 

C  READ  FIXED  VOLTAGE  CARDS 
NBDY-0 
NTCT-NX*NY 
.NEL- ( NX- 1 )  •  ( NY-1 )  *2 
NBK-0 

DO  30  II-1,NEL 

I?(II.LS..NTOT)  W(ii)-o. 

C-iD(  II )  =0  . 

80  CONTINUE 

DO  so  II*1,NV 

READ ( 4. «, END- 100)  VS, II , 12 , J1 , J2 , ILB 
ISdl.EQ.l)  VF-VINN 
IS(II.SQ.2)  VF-VCUT 
IS(IL3.EQ.9)  THEN 
NBK-NBK+1 
I3K(NBK)-I1 
SNDI? 

DO  60  JJ-J1,J2 
DO  70  XK-I1,I2 
.NBDY-NBDY+1 
NN-NX*(GJ-1)+KK 
MBDY(1,NBDY)-MN 
MBDY(2,NBDY)— ILB 
IF(ILB.EQ.9)  THEN 

VP ( NBDY ) - VInN-SXON/RKON* ( 1 . /RIP- I . /RR ( NN ) ) 
VV(NN)-VP(NBDY) 

ELSE 

VP(NBDY)-VF 

VV(NN)-VF 

ENDIF 

70  CONTINUE 

60  CONTINUE 

50  CONTINUE 
RETURN 

100  CONTI.NUE 

STOP  'NXTIN' 

END 

SUBROUTINE  VEXP ( NRP , IN , RR ) 

IMPLICIT  REAL'S  (A-H.0-2) 

DIMENSION  IN(*),RR(*) 
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n 


'c  :3  = 

?.R  '  II  '  )  -RR(  II  ^ 

IR-  RR;  II)-RR(II-I';  ) /FLOAT!  IN(  II  ;-IN(i:-l)  ) 
ic  10  ij=»:n(II',-i.im(ii-i;*i,-i 

?.R  ,  J  J  '  »  RR  (  J  J I  I  -  DR 
10  ic.'rri.'ruE 
10  CCNTI.'RIE 
RETURN 
ENC 

EUBRC'JTINE  SLIPS  (  EGI ,  RRI ,  3EI ,  THI ,  DEL  ,  EGO ,  RRO  ,  3E0 ,  THO ,  ITEST ) 
IMPLIEIT  REAL»9  (A-H,0-Z^ 

LCGIC.\L  3MX  ,  QMN  ,  ^TP  ,  QBT ,  QTC ,  QBC ,  QEX 
CCMMCN  /INST/  RI?,RCP,RIA.ROA, SNOT, EG2 
SAVE 

DATA  NCALL  /O/ 

IF'NCALL .SQ. 0 )  THEN 
PI-ACOS.'  -1  .  ) 

?i:2»: . 5*?I 
RTCD-130. /PI 
;"'ALL-1 

RKCN- ( 1 . /RIP- I . /ROP ) 

SKCN-ENCT* ( ROP+RI? ) ‘RKON 
VINN*SKCN/{1 .-RIP/ROP) 

VOUT^vinN-EXCN 
liQ^SKLN/RKCN 
EKuIF 

TKE-90 . -VHI 
RDEL=DEL/RTCD 
TSET  =  TAN  ( 3EI/RTOD ) 

RTKE»rKS/RTOD 
EKC»CQ/(2.*RRI) 

CC=SXC/EGI/CCS(BEI/RTOD)**2 
•J'JRT=CC -T3ET*S1N  f  RJDEL )  +  ( 1 .  -CC )  *COS  ( RDEL ) 

RR  =  RRI/UURT 
?MAX»?I/2. 

:?!CC.NE.:. )  ?MAX»-ATAN(TBET/(1.-CC)  ) 

'JUPM«CC-T3ET*SIN  (  PMAX )  +  ( 1 . -CC)  »COS  ( PMAX ) 
i;0’CM-CC-T3ET»SIN  (  PMAX )  -  ( 1 .  -CC )  *COS  (  PMAX ) 

U1 8  0 -CC-T3ET-*SIN  ( PI-RTHE )  +  ( 1 . -CC )  *COS  ( PI-RTHE ) 
U270-CC-T3ET*SIN ( PI 32-RTHE )  +  ( 1 . -CC) *COS ( PI32-RTHE ) 
RMAX-RRI/UUPM 

rmin-rri/uucm 
R180»RRI/Uia0 
R270=*RRI/U270 
PMAX*  ?MAX*RTOD+THE 
Ir!?MAX.LT.O. )  PMAX-PMAX+360 . 

I?(?MAX.GT.36C. )  ?MAX»PMAX-360. 

ITEST  LOGIC 

IFfRMAX.GT.RMIN)  THEN 
PMIN-PMAX+iaO . 

ELSE 

RT*RMAX 
RMAX-RMIN 
?T  '^MAX 
PMAX-PMAX+130 . 

RMIN=RT 
?MIN=?T 
SNDIS 

I?;?MJ\X.LT.O.  )  PMAX-PMAX+360. 
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IF  :  PMIN.IT  .  0  .  )  . 

IFfFMAX.UT.  360  .  1  ?MAX»?HA.X-  3  60  . 

IF'  PMI.N.3T  .  360  .  )  ?jlI.M-?MIN-3  60  . 

:;mx=  ?kax  .  gt  .  the  .  and  .  ?max  .  ge  .  the^del  .  and  .  fmax  .  ge  .  ecp 

QMN=  -=MIN .  GT  .  THE  .  AND .  ?MIN .  GE  .  TI-E-^DEL  .  AND  .  EMIN  .  IE  .  El? 

if'Qmx.and.gmni  then 

GMX=CMX . AND . ?MAX . IT . EMIN 
';mN=  .  NOT  .  QMX 
END  I? 

QT?=QMX . OR . ERI . GT . EC? . OR . ER . GE . SOP 
QBT=CHN . OR . ERI . GT . RIP . OR . ER . GE . El? 

CTC  =  THE.GT. 130  - .AND.THE+DEG.GE. 130 . . AND . El 80 . GE . ROA . AND . R1 30 . GT . 
tROP 

gBC  =  THE.GT.  lao  .  .AND.THE-t-DEG.GE.  130  .  .AND.R130  .GE.RIA.AND.R130  .GT. 
-RIP 

CEX  =  THE.GT.27  0  .  .  AND  .  THE-t-DEG  .  GE  .  27  0  .  AND  .  R27  0  .  GT  .  RCP  .  AND .  S270  .  GT .  RIP 
C 

2TKE»RDEG 

TT1=MIN(MAX(0 . ,RDEL) ,?I-STHE) 

TT2«MIN(MAX(PI-RTHE,DEL)  ,Pr3  2-R'rHE) 

IF(CT?.OR.QTC)  THEN 

C  HITS  TCP  PLATS  WITHEN  THE  +  DEL  INTERVAL 

IF(CTP.AND. (PMAX.LE.130. . OR .S130 . GT .ROA ) )  'THEN 
ITEST«=2 

CALL  NEWT(CC.T3ET,TT1,RRI.R0P,2THE) 

ELSE  I?(0TC)  THEN 
ITEST-4 
2THE»?I-RTHE 
ELSE 

ITEST=»6 

CALL  .NEWT{CC.TBET,TT2,RRI,R0P,2THE) 

ENDIF 

ELSE  IFIQBT.OR.QBC)  THEN 
C  HITS  BOTTOM  PLATS  WITHEN  INTERVAL 

IStQBT.AND. ( PMIN. LE. 180. .OR. R180.lt. RIA) )  THEN 
ITSST-3 

CALL  NEWT (CC.TBET.TTl.RRI, RIP, 2THS) 

ELSE  IF(QBC)  THEN 
ITSST-S 
2THE-PI-RTHE 
ELSE 

ITSST-7 

CALL  .NEWT  CCC,T3ET,TT2,RRI,  RIP,  2THE) 

ENDIF 

ELSE 

ITSST-l 

ENDIF 

UURT»CC-TBET*SIN ( 2THE ) + ( 1 . -CC ) *COS { 2THE ) 

3EIN- ( 1 . -CC ) *SIN ( 2THE ) +T3ET*C0S ( 2THE ) 

BEO-0 . 

IF '  HURT  .  NE  .  0  .  >  BEQ«.\T.AN  f  3EIN/UURT 1  *RTGD 

V  j  J  O  O  . 

IF(3EO.LT.-90. )  BEO-3EC+360. 

RRO-RRI/UURT 
THO-THI- 2THE*RTOD 

SGO-EG2+VINN-S.KON/RKCN*  ( 1  .  /RIP- 1 .  /RRO ) 

RETURN 

END 

SUBROUTINE  FWKM ( NE , lESM , RXST , EKINC , SBAR , S IG , ERAN ) 

IMPLICIT  RFAL*8  (A-H,C-Z) 
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max:<=»:esm(  i ; 

EETERMI.'rE  ZS.\3. .  SIGI-IA 
.'rN=  j 
E3AR=0 . 

■'SGM-  D 
SQSM«0  . 

EC  13  ::»i..'rE 
nn»nn*:esm( ;r 1 

YV-SKST-^i  :r-l  1  *EKrMC 

ysum-;si:m-i-''y»iesmi  id 

SQSM-5CSM+YY*';''*:zsM( ::  > 

10  cc.'Trr.vuE 

lE'NN.CE. 3 )  THEN 
EBAi?-'S'JM/NN 

5IG- ( SCSM-EaAR*SBAH*NN) / (NN-1 ) 

:e;s:g.gt.o.  )  sig-scrtoig) 

END  I? 

MAXCH-1 

ha:<x«:ssm(  : ) 

DO  so  :»2.NE 

IE( lESMf T ; .GT.MAXX)  THEN 
max;<>;esm(I) 

MAXCi-I 

ENDIE 

so  CONTINUE 

HMAX « ELCAT { MAXX ) / 2 . 

CHl-0. 

CH2-NE 
DEN»1 . 

DO  SO  :«1,MAXCH 

:?!HMAX.GT.?LOAT(IESM(I) ) . AND. HMAX.DS. FLOAT ( lESM ( 1+ 1 ) ) 
DEN-IESM( I+l ) -ISSM( I ) 

IF(ABS(DEN) .LT.l .E'S)  DEN-1. 
CHl>:-l+(HMAX-ELOAT(IESM(I) ) )/DEN 
ENDI? 

60  CCNTI.NUE 

DO  70  I-NE,.MAX(2,.MAXCH)  ,-l 

IE;HMAX.GT.?L0AT(IESM(I) ) . and. HMAX.LE. float (IESM(I-I) ) 
den-;esm(1-i)-iesm(I) 

CH2-:-l-(HMAX-?LOAT(IESM(I) ) ) /DEN 
ENDIF 

70  CONTINUE 

CCH1-EXST^CH1*EKINC 

CCH2-BKST+CH2*SKINC 

ERAN- ( CH2-CH1 ) *EKINC 

RETURN 

END 

SUBROUTINE  NEWT (CC.T3ET. DEL, RRI.RRT. THE) 

IMPLICIT  REAL'S  (A-H,0-Z) 

UU ( TT ) -CC-TBET*SIN ( TT ) + ( 1 . -CC ) 'COS ( TT ) -URRT 
UP;TT)--TBET-»C0S(TT)-(1.-CC)'SIN(TT) 

IFRTY-0 
GSS-4 . 

40  CONTINUE 

URRT-RRI/RRT 
X?P-DEL,'GSS 
E?P-UUf X?P) 


)  THEN 


)  THEN 
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;<?^x??-;p?  OYPP 
DO  -.3  :r>:.;o 

=  I  X? ' 

OYP^'JP'XP' 

OD»xp-:<pp 

AA=  ; . "OYP^Ovpp- 3 . * , v?-v?p  < /CD  I / OD 

3B=vp , 3yp 

tH£=;<P-3B*  '  1  .  -BB».\A/DY?1 
iRS-ABS' . THE-XP'/THS; 

II  .  3T  .  :0  .OR  .  lERTY  .2Q.  I  PRINT*  , 
lE'ERR.lT. - .E-i )  30  TO  20 

;<pp=xp 


ERR, THE  -  ',11, ERR. THE 


:<P=THE 
vpp,  /p 

0YPP-3V? 
10  CONTINUE 
30  CONTINUE 


3SS-GSS»2. 

C  PRINT* II , ERR, THE  -  •  .  II , ERR , THE , GSS 
I?'IERTY.CS.2)  GO  TO  40 
RETURN 

20  CONTINUE 

I? ' THE . GT . OEL )  GO  TO  30 

RETURN 

END 

SUBROUTINE  SETXO  ( NODE  .NBAND , NEL .  I7AXIS  ,  MOON ,  SX  ,  XX ,  YY  ,  ANEL ) 
IMPLICIT  REAL'S  (A-H,0-Z) 

C 

C  THIS  SUBROUTINE  SETS  UP  MATRIX  XI (Jl) 

DIMENSION  XX(») ,YY(*) ,SX(NBANO,*) ,ANEL(*) 

DIMENSION  MC0N(7,*) 

DIMENSION  3AR(0;3) .3(3) ,C( 3) ,XE( 3) ,YS( 3  ) 

C 

C  CLEAR  SX 

DO  10  I»l,NODE 
DC  10  J-1, NBAND 
SXf J,I)-0. 

10  CONTINUE 

C 

C  ASSEMBLE  K  INTO  SK  IN  PACKED  FORM.  RHS-Q*CHD  INTO  W 
C  LOOP  OVER  ALL  ELEMENTS 
DO  200  IE- L, NEL 
3AR( 3)-0. 

DO  210  1-1,3 

XE(I)-XX(MCON(I.IE) ) 

YE(I)-YY(MCON(I,IE) ) 

BAR( 3)-BAR(3)+XE(I)*XE(I)+YS(I)*YE(I) 

210  CONTINUE 

3(1)-YS(2)-YE(3) 

3(2)-YE(3)-YE(l) 

3(3)-YE(l)-YE(2) 

C(1)-XE(3)-XE(2) 

C;2)-XE(1)-XE(3) 

C( 3 )-XE(2)-XE(l ) 

C 

C  LOOP  OVER  NODES  OF  ELEMENT 
BAR(0)-1. 

3AR(l)-(XE(l)-fXE(2)+XE(3)  )/3. 
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on  no 


3AK  ; i -££( 3 ) : / : . 

3AR:  3 ; -SQRT' 3AR! 3 ' . 

A-;»4  .  '.vrcii  /  BAR;  ITAXIS'I 

3C  ;so  :"i,3 

:N=MCCNf I , lil 

30  :so  J-: . 3 

3N-MCCN(J.;S; 

IN. ON! 

I'.N-IN-IN+I 

IE' JN.OT.INl  :<n-on-:n+i 
ElEX-  3II;*3( Ji*C;Ii»C;J) )/A4 
3X  :<N,0N)=3X(;<N,0N)-£LZX 
250  CONTINUE 
200  CONTINUE 
RETURN 
END 

SUBROUTINE  SETQ 2 ( NODE , NEL . DE3EE , Q , G , ANEL , MCON ) 
IMPOICIT  REAL»3  (A-H.0-2) 

C  THIS  SUBROUTINE  SETS  UP  THE  RHS  VECTOR  Q 
DIMENSION  Q(*} ,0(*1 .ANEL( » ) ,MCON( 7  ,  » ) 

CCZAR  3 

DO  10  I-1..NCDE 
G(I)-0. 

LQ  CONTINUE 
DS3I-1 . 

I?(0E3iS.CT. 0 . )  DSQI-I ,/DEBYE/DEBES 

COMPUTE  NODAL  AREA  AND  CHARGE 
DO  20  lE-l.NEL 

A3-ANEL(IE)/3.*DSQI 
DO  20  1-1,3 

:n-mcon(I,is) 

Q(IN)-g(IN)+G(IE)*-A3 
20  CONTINUE 
RETURN 
END 

SUBROUTINE  SSSBC ( SX . Q , ED , NRD , NODERD , NODE , NBAND ) 
IMPLICIT  REAL»9  (A-H,0-Z) 

DIMENSION  SK(NBAND,*) ,Q( • ) ,ED( * ) ,NCDERD( 2 , * ) 

DO  10  I-l,NODE 
DO  10  IRD-1,NRD 
:<- NODERD  (l.IRD) 

IE(AaS(:<-I)  .LT. NBAND)  THEN 

IF!X.GT.I)  Q(I)-Q(I)-SK(X-I-t-l,I)*FD(IRD) 
I? ( X . LT . I )  Q (I ) »Q ( I ) -SK ( I-X+1 , K ) *FD ( IRD ) 
END  I? 

10  CONTINUE 

DO  20  IRD- 1,  NRD 
K-NODERD(l.IRD) 

SK( 1 ,X)-1. 

Q(X)-?D(IRD) 

DO  30  J-2, NBAND 
SX(J,K)-0. 

IE(X-J.GE,0)  SX(J,X-J+l)-0. 

30  CONTINUE 
20  CONTINUE 
RETURN 
END 

SUBROUTINE  SOL'/ '  A  ,  3  ,  NEQT , NBAND  ) 
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IMPLICIT  aEAL’8  (A-H,0-Z', 
DIMENSION  AiNBAND.*) ,3(*) 
NSCLVE-NEGT-1 
NB-NBAND-I 
J;»NECT-NB 
DC  IQO  I>:,NSQLVE 
G=i.o/A(i,i: 
lE'I.OT.Cr;  NB»NEQT-I 
NB1»NB+I 
DO  200  II»:,N3 


s-A(ii-:,r)«Q 

I?'ABS/R) .3Q.0. )  GO  TO  200 
3(  JI)-8(CI)-R*3(I) 

”»nb:-ii 
DO  300 

AdJ.JD-Adu.JIj-R^Adl-f 
300  CONTINUE 

200  CONTINUE 

100  CONTINUE 

NB-NBAND-1 
DO  400  I-1,NEGT 
Q=l./Add) 

3d)-3d)*Q 

lE(I.SQ.NEQT)  GO  TO  400 
IB-MIN(NB,NECT-I) 

DO  500  ri-l.IB 

AdI+ld)-A(II-dd)*Q 
500  CONTINUE 

400  CONTINUE 

DO  500  r-l,NSOLVE 
J-NECT-I 
IB«MIN(NB.I) 

DO  600 

3(  J)«B(J)-A(II+l,J)*3(II-t-J) 
600  CONTINUE 
RETURN 
END 
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on  o  o  o  o 


?P.CCKMi  ;-iEMIV 
iMFiicrr  REAi:,*3  ;a-h.o-:) 

Cr-IAFACTEH  TESTt  101 , 131  I  .ICHARdO)  ,1TMP'7)*14 
IIME.VSION  :TCT;  101  .  lOl  l  ,.'fSUM(  101  1  dSCJMf  101  !  ,JSUM(  101  )  , 

-P.HCEf  15  I  ,13’JM(  101  ,  .MS'JM(  101 1  ,NECT;  101 1  ,NAOT(  101 ) 

3ATA  rCHAR/ ' 

3ATA  ISu'M. JSUM.LSUM,MS'OM,.'JSUM,NECT,.SAOT  ,707*0/ 
lATA  TEST  /102Q1*' 

3c  10  :=i,4i 
OC  10  J=l,41 

:tqt;i,J)*o 

cc.'rriNUE 
?l=ACOS( -1 . ) 

3An-?I/130 . 

TKE'rA»90  . 

THETA*  THETA*  SAD 

st»3:m(Theta; 

M13Q=0 

SET  QUADRATURE  LIMITS 
ME  =  41 
MP  =  4l 
MR- 11 
MA-11 
ISW6-1 

I? ( :SW6 .EQ. 1 )  THEN  _ 

PRINT*, 'M130,ME,MP,MR,MA  -  *  , M130 ,ME ,MP ,MR ,MA , '  ENTER  NEW  VALUE 

-S' 

READ(5,*)  M130,ME,MP,MR.MA 

PRINT* , 'M130 ,ME,MP,MR,MA  -  ' ,M130 ,ME ,MP ,MR ,MA 
ENDI? 

DEFINITION  OP  PLATS  AND  APERTURE  RADII 

DEFINE  GEOMETRY 
ENOT-1000. 

CALL  GECMA ( VOUT , VINN . SIP , SOP . SIA , SOA . SNOT ) 

DEFINITION  OF  PLATS  AND  APERTURE  RADII 
RB»(RI?+ROP)/2 
DRP-ROP-Rl? 

DIA-ROA-RIA 
RA-DIA/2. 

EAL-DIA 

SXCN-ROP*RI?/ ( ROP*RCP-aiP*RIP ) 

VI- 1 , / ( EXON* ( 1 . +RIP/ROP ) ) 

VC-VI-1 ./EXON 
QO- ( VI-VO ) / ( 1 . /RIP-1 . /HOP ) 

SET  130  DEG  APERTURE  PARAMETERS 
IF(M180.E(3,1)  THEN 
C  SLIT  OF  CONSTANT  DELTA  R 
RC-EAL/2. 

WRITE ( 6, ' (lOX, ’ 'APERTURE  AT  180  IS  SLIT.  WIDTH  -'',F8.3, 

-  '  '  CM '  '  )  '  )  EAL 

ELSE  IF(M180.EQ.2)  THEN 
C  CIRCLE  OP  RADIUS  RC 
RC-RA 

WRITS ( 6, ' (A1,6X, ' 'APERTURE  AT  180  IS  CISCLS.  RADIUS  -'’,F8.3, 
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F 


h:i3z 

ZZZAL'LT  ;S  SO  APERTURE 
M130-0 

RC= I RCP-RI?'/0. 

>fRITE,6,'  SX.''MO  APERTURE  AT  130  OEGREES  ' 
ESDI? 


C 

C 


C 

c 


c 

c 


OEErSITIOS  OF  U  QUADRATURE  P.AHAMETERS 
UI?=Ra/RI? 

UOP-SB/ROP 

JIA»aB/RZA 

UCA=.RB/ROA 

UINC*  j  UIA-UOA ) /MR 

UST»UOA+UISC/2 

:f(mr.ie.1)  then 

MR=1 


UST=1. 


ENDI? 


DEFINITIONS  OF  ALPHA  QUADRATURE  PARAMETERS 
ALF3 1 - - AS IN ( RA/RB ) 

ALFB2»-ALFBl 

ALFINC- ( ALF32-ALF31 ) /MA 

ALFAST-AL?Bl*AL?INC/2 . 

I?(MA.LE.l)  THEN 
MA*1 

AL?AST«0 . 

ENDI? 


DEFINITION  OF  EX  QUADRATURE  PARAMETERS 
SK1-0.3S 
SX2=1.15 

SKINC= ( EX2-SK1 ) /ME 
SFUD-.0099 
SFUD«0 . 

EKST-EKl+EKINC/2 . +EFUD 
DE»(VI-VO)/21. 

EB-1.-10.*DE 
IF(ME.LS.l)  THEN 
ME»1 
EKl-1. 

EX2=1. 

EKINC»1. 

EKST-l.+EFUD 

SNDIF 


DEFINITION  OF  PHI  QUADRATURE  PARAMETERS 
PHll  — 12.0/ST 
PHI2-12.0/ST 
PHI1-PHI1*RAD 
PHI2-PHI2»RAO 
PHINC - ( PHI 2 - PHI 1 ) /MP 
PHI  ST»  PHI  1 -f  PHINC/ 2 
I?(MP.LS.i)  THEN 
MP-1 
PHIl-0. 

PHI2-0. 

PHINC- 1. /ST 
PHIST-0. 

ENDI? 
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AS'JM-  : 

RIA-S3-?.C 

HCA-R3*3C 

•jia»s3/h:a 

•JCA-S3/P.CA 


r  2G  ^UADBA7'Jli£S 

y!MR»  MK-  -  /  2  -  : 

MMA-  HA-. 

MAT=: 

iF'Miao ..'rs. ; '  then 

MAT«MA 

MMA-: 

MA-: 

MB-MAXfMA.MAT ; 

.'JTSJ’-O 
30  30 

U-CJST-  I  3-1 )  «aiNC 
3*a3/U 

rE'MR.EG.l)  •JIMC»U*U/RB 
;?(MB.3S.l)  AiFINC-U/RB 
C 

C  DEFINITIONS  OF  ALPHA  QUADRATURE  PARAMETERS 
0R»A3S(RB-R) 

CALF- I R*R+R3*Ra-RA*RA) /2 . /RB/R 
ALFA1--AC0S(MIN(1 . ,CALF) ) 

ALFA2--ALFA1 
ALFA-0 . 

C 

DO  30 

IF'Miao .EQ.2 )  THEN 

ALFA- ALFAST+ ( MM- 1 ) *AL?INC 

I?(ALFA,LT.ALFA1.0R.ALFA.GT.ALFA2)  GO  TO  30 
SQ-SQRTCRC^RC-CRB^SINCALFA) )**2J 
RCAL-R3*CCS ( ALFA) 

RIA-RCAL-3Q 

ROA-RCAL-fSQ 

UlA-RB/RIA 

•JOA-RB/RCA 

ENDIF 

C 

DO  40  I-l.ME 

SX-EXST+ ( I-l ) *SKINC 
SGT-SXtVI-QQ* ( 1 . /RIP-1 . /R ) 

DO  40  I-1,MP 

?HI-?HIST+( J-l)*PHINC 

C  CALL  TRAJEC(Ra,U, ALFA, THETA, PHI, EX, ITEST,3H0,UIP,U0P,UIA,U 

C  +  OA) 

NTRJ-NTRJ’+l 

CALL  SLL?S(EGT,R,  PHI/RAD,  0.00,270.  DO,  QQ,  RIP,  R0P,R1A,R0A, 

.r  RR,  BET,  ?MAX,RMAX,3MIN,R180,  REND,  ITEST,NTRJ) 

SSND-EX-rVI-QQ*  ( i .  /Rlp-1 .  /RR ) 

DO  SO  M-L,.MAT 

IF(M180 .NE.2)  THEN 

ALFA-ALFAST+ ( M-1 ) *ALFINC 
E.NDI? 

IF(MR.EQ.l.AND.MB.SQ.l)  THEN 

•  1  ^  T 
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ME  .  .'XE  .  1  ;  THE.V 


ZE'MP.NE.:  THEN 

1 2=«J 

ELSE  :E'MR..'IE.1  ‘  THEN 
I  L»  L 
ELSE 

;l=«ma;<  '  M  ,;-;m  ] 

ENDI? 

ELSE  r??MF.:rE.:!  then 
ZL=T 

I?'MR..'JE.  Z  :  THEN 
IZ»L 
ELSE 

LZ»MAX(M,MM) 

ENDIE 

ELSE 


L2-MAX(M,MM) 

MMR»MR 

ENDI? 

TEST(ri.:2).'Q* 

IF(M130. NE. 2. AND. (AX.EA.lt. ALEAl .OR. ALFA. GT.ALFA2) )  GO  TO 
+  60 

TEST  (  ;i ,  r 2  )  -ICiAR (  XTEST ) 
lEdTSST.NE.i;  GO  TO  SO 
ISUM(I)-ISUM(L)+: 

JSUM  ( J  ) -SSUM  ( J  )  + 1 
LS'JM(L)-LSUM(L)*1 
MS'JM(M)-MSUM(M)-: 

ASUM-ASUM+1 . /U/U/U 

:tot(I,J)«itot(I,J)+i 

.NRHO-  ( REND-SI? )  /  ( ROP-RI? )  *21 .  +1 . 

.VRHC-MIK(MAX(  I  ,NRHO )  ,  21 ) 

NSUM ( NRHO ) -NSUM ( NRHO )  + 1 
NEO-(SSND-VO-l. )/(Vl-VO)*21.+l. 

.NEO-MIN  ( MAX  ( 1 .  NEO ) ,  2 1 ) 

.NAO-(BET+8.  )/16.*21.t-l. 

NAO-MIN(MAX(l,NAO) ,21) 

NEOT ( NEO ) -NEOT ( NEO )  +1 
NAOT ( NAO ) - NAOT ( NAO ) + 1 
60  CONTINUE 

40  CONTINUE 

1P(L.EQ,MMR.AND.MM.EQ..MMA)  THEN 
IF(Il.GT.l.AND.I2.GT.l)  THEN 
DO  120  JJ-L,.MIN{  41,12) 

WRITS(6, ' (1X,41A3) ' )  (TSST( I , JJ ) , 1-1 ,MIN( 41 , II ) ) 


120 

CONTINUE 

ENDI? 

ENDI? 

30 

CONTINUE 

IF( 11 .GT . 1 . AND . 12 .GT . 1 ) 
WRITE (6, 33)  PHIl/RAD 

THEN 

33 

FORMAT (//•  PHI  START- 
DO  80  II-1,MIN(41.MP) 

'  ,FS.l, IX, 'DEGREES' ) 

WRITS(6, • (IX, 4113) ' 

)  (ITOT(I,II) ,I-1,MIN(41,ME) ) 

80 

CONTINUE 

WRITS (6, 34)  PHI2/RAD 

34 

FORMAT (IX, 'PHI  END-', 
WRITE(6,3S)  SK1,SX2 

F4.1,1X, 'DEGRESS' ) 

102 


35 


35X, 'ENERGY  END-',??. 5) 


rCRMAT  •  :x  .  ■  E.‘:ERGY  start-  ■  ,  E ’  .  =  , 

E-'ro:? 

'■^ETE  .'5,1001 

I'OC  rCRMAT:  /'  9X,  '  E3EG  '  .  ll.X,  'PHI'  .14X,  'R'  .IIX,  '  .U.?  '  ,liX,  'SEND*  ,11X, 
-  3E.'rE  '  ,  i;x  .  '  ?MC  ' 

ZC  14  :  =  :  .max; MR .MB .ME.MP, ;i ) 

ITMP ' 1 ; » '  ’ 

ITMP ' 1 ; - '  ' 

ETMP'3l-'  ' 

ITMP '41-' 

ITMP ' 5 . - ' 

ETMP ( 6 ) » '  ' 

ITMP' 7 '  ' 

e:<-exst-  ( r-1 )  'exinc 
?KI-?HI3T*( :-l ; *?HENC 
3PHE-PHI/RAD 
y-usT-(r-i)*uiNC 
R-1  ./U 

AEEA-ALE.XST*  ( I-l )  *AI,?INC 
EC-E3-;-?ECAT(  I-l )  *DE 
3C--3.-(ELCAT(r)-.5)»16./21. 

RC-Rr?+ ( ELCAT ( I ) - . 3 ) • ( RCP-RI? ) /21 . 

r SS-  IS’JM  ( I )  +JSUM  ( I )  +ESUM ( I )  +MSUM  (  2  )  +NECT  ( I )  +NAOT  ( I )  +NSUM ( I ) 

:e(iss.ie.O)  go  to  24 

:e(I.es.me)  writs (Ltmp( 1) ,4001 )  isuM(r),sx 

lEd.EE.MP)  WRITE(LTHP(2)  .4001)  JSUM(I),3PHI 
lEd.LE.MR)  WRITS(LTMP(  3)  ,4001)  LSUM(I),R 
I?(I.EE.MB)  WRITE(LTMP{4)  ,4001)  MSUM  ( I )  ,  ALFA/'RAD 
IE(I.2E.21)  WRITS(LTMP(S) ,4001 )  NECT(I),EO 
:?(r.EE.21)  WRITS(r-TMP(6)  ,4001)  NA0T{I),30 
:E(I.2E.21)  WRITS(ETMP(7) ,4001)  NSUM(I),30 
4001  FORMAT (14, no.  5) 

WRITS (6, 3)  LTMP 
24  CONTINUE 

5  ?CRMATC7(1X.A14) ) 

C  ’WRITS(6,101) 

DO  23  1-1,11 

IISM-:iSM+ISUM(I) 

23  CONTINUE 

GF-  (  ASUM )  »R3*RB»ST'*ST*PHINC'*UINC<EKINC'*AL?INC 
THETA- THETA/RAD 
'WRITE(6 .102) 

102  F0RMAT(1X/1X/1X, '***TKETA*-*SUM  OVER  AL’  GSCMETRIC  F.aCTOR '  ) 

•WRITS  (5,9)  THETA ,  IISM  ,  GF 
9  FCRMAT(1X,F10.5,I10.1PS1S.5) 

IF(ME.GT.2)  THEN 

CALL  FWHM(ME,ISUM.EKST,SXINC,SBAR.SIGE,SRAN,MAXX) 
ERAN2-FL0AT ( IISM ) •EKINC/FLOAT ( MAXX ) 

2001  FORMAT ( //IX, 'SBAB, SIGH,  DELTA  S  OVER  E  -  ',4F10.5) 

WRITS (6, 2001)  E8AR,SIGE.ERAN/MAX(EBAR,1.D-S) ,ERAM2 

SNDIF 

IF(MP.GT.2)  THEN 

CALL  FWKM ( MP , JSUM , PHIST , PHINC , P3AR , SIG? , ?RAN , MAXX ) 

PRAN2- FLOAT ( IISM ) »PHINC/FLOAT (MAXX) /RAD 
'WRITS ( 5 , 2002 )  ?3AR/RAD,SIG?/RAD,?RAN/RAD,?RAN2 

2002  FCRMAT(//1X, ’?3AR,SIG?,PFWHM  -  ’,4F10.5) 

SNDIF 

IF(.MR.GT.2)  THEN 

CALL  FWHM  ( MR  ,  LSUM ,  UST , UINC ,  UBAR ,  SIGR ,  URAN , .MAXX ) 
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■jran;-?:.oat;  :rsM)  *'jinc/ float  .'ma>cx  . 

>m:TF : 6  .  0Q0  3  1  : .  /  UBAH .oXGS , ORAN , 0RAN2 
;0C3  FORMAT'/ 'RBAR,  3  IGR.RFWHM  •  .iFlO.S) 

F.'fCFF 

IF  MB. GT. 21  THEN 

G-ALL  FWKH  (  MB  .  .MSGM  .  ALE  AST  .  ALE  I.NC .  .ABAR  ,  3  IGA  ,  ARAN  ,  .'iAXX 
.ARAN2  »  FLOAT  '  I  ISM  1  •  ALEINC/  FLOAT  <  MAXX /RAO 
•VRITE.  5,2004  1  .ABAR/ RAD  .  5IGA/ RAD  ,  .ARAiN/RA  0  ,  ARAN2 
2CC4  FORMAT;/  /:x.  .ABAR.SIGA.AEWHM  «  ',4F10.5) 

END  IE 

->  e  -  _  1  ^ '  2 1 

E3« . 34-OE/ 2 . 

0.VL1  r/mM  (  21  ,  .NECT  ,  EB  ,  OE  ,  SCBR  .  SIGE  ,  EORAN  ,  -MAXX 
•/^ITE;  6  ,  2005  )  EOBR,  SIGE,  SCRAN 

2005  FORMAT'//IX, 'EOBR,SIGE,SFWKM  »  ',3F10.3) 

03»:5  ./21  . 

3B=-3 . -DB/2. 

CALL  r//HM(  21  ,.VACT.Ba,DB,3CBR,SIG3,3CRAN..MAXX) 

WRITS! 6 , 2006 )  BOBR , 3IG3 . 3CRAN 

2006  ECRMAT(//1X, *3CBR.SIG3,3FWHM  -  ',3F10.5) 

DR-(RCP-RI?;/21. 

RHCB-RI?*0R/2. 

CALL  FWHM{  21  ,.NSUM.RHCB,DR,ROBR,SIGR,RCRAN,.MAXX) 

WRITS.- 6 , 2007  )  ROBR.SIGR.RCRAN 

2007  FORMAT (// IX, 'R0BR,SIGR,RFWKM  -  ',3F10.5) 

END 

SL'BRCL’TINE  TRAJEC  ( R3  ,  G  ,  ALEA  ,  THETA ,  PHI ,  EK  .  ITEST ,  RHO ,  DIP  ,  UOP  ,  UIA ,  UOA 
-) 

IMPLICIT  REAL*3  (A-H.C  Z) 

SAVE 

DATA  .VCALL/0/ 

I?{NCALL.NE.O)  GO  TO  10 
PI«ACCS(-1. ) 

TPI*2 . *?I 
^«?i/13Q. 

RAD270-270.*RAD 
C3-C0S ( 0 . *RAD ) 

RHC-0. 

S3-3IN(0.*RAD) 

10  CCNTI.NUE 

.SCALL-NCALL+1 
W-SX+2* ( U/UIP/UOP-1 ) 

C  CALCULATE  TRIG  QUANTITIES 
CA-3IN(ALFA) 

CA-COS(ALFA) 

CT-COS( THETA) 

ST-S1N( THETA) 

SF-31N(PH1) 

CF-COS ( PHI ) 

C  MAKE  ANGULAR  TRANSFORMATION  TO  THE  TRAJECTORY  PLANE 
XI— 3A*ST»SF+CA*CT 
Y1-3T»C? 

G-SCRT(X1*X1+Y1»Y1 ) 

SG-Xl/G 

CG-Yl/G 

X2-CA*ST*SF+SA»C7 

Y2--3G*CA»CT+SG*SA*ST*SF-ST*CF*CG 

C2W-Y2*Y2/(X2*X2+Y2»Y2) 

C  CALCULATE  ORBIT  CONSTANT 
UP»U*U/W/C2W 
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O  fl  o 


is;<=u»x2,  y; 
icx-i;- jp 

-xr?-ATAN: ; hsx . £cx 
:?5:?S»3GRT  :  ESX»3SX-£CX«SCX'. 

:£=■ 

:<A»  'JCP-JP' /'JPEPS 

■  AflS  ;  :<A  ;  .  E£.  :  .  THEM 

:c;^»Accs.  ;<A , 

rHiA=:tA-c3i? 

e:-;i3»-:<a*chi? 

ZE'CHIA.lt.O.  E:-::A=~-:IA-rT?: 

;e'CHT3.et.  : .  c:xi3-c:-:i3-^t?i 

:?'C-;iA.3T.rp: )  EHrA«cHiA-T?r 
:E'C-iI3.3T.T?:  :  CH:3-CHI3-T?I 

IE ' c:-:r 3 . 3e .  . ame  .  chib  . it  . rad2to  )  cHiA-c-as 

3I3TA* 0  . 

EMEIE 

X3-  'JTP-UP  )  /UPEPS 
:E(AaS(X3) .32.1. )  THEM 
X3-ACCS{Xa) 

CXIC-XB+Cil? 

c:-ii3--xB*c-a? 

:?;aic.iT.o. )  c-iic*CHiC4-TPr 
:f{C-ii3.3t.o.  )  chib-ctio+tp: 

:E:CTic.GT.TPr chic»chic-tpt 
;?;CTi3.crT.TPi)  ctio-chid-tpi 
:?(C-iI0.3E.?I  .ANC.CHI3.3T.RAD27  0)  CHIC-CHID 
2332 

CHIC-0 . 

2MDI? 


ICES  THE  CR3IT  GO  THBCUGH  THE  AMALI2ER? 

IS  THE  FIRST  EXTREMUM  A  MAXIMUM  OR  A  MINIMUM? 

IF'C3I?.3T.O. .AMD.UIP.3S.UP+UPEPS)  THEN 
C  A  MAXIMUM  -  DOES  IT  HIT  THE  INNER  PLATE? 

ITEST-3 

C  A  MINIMUM  -  DOES  IT  HIT  THE  OUTER  PLATS? 

ELSE  IFfCHIP.LE. 0. . AND. UP-UPEPS . 32. UOP )THEN 
ITEST-2 

ELSE  lE'US.CT.UI?)  THEN 
ITEST-3 

ELSE  IF(UE.3E.UOP)  THEN 
ITEST-2 

C  DOES  IT  GO  THROUGH  THE  EXIT  APERTURE  OF  THE  HEMISPHERE? 

ELSE  IF(UE.GT.UIA)  THEN 
ITSST-3 

ELSE  I?(UE.3S.UOA)  THEN 
ITEST-4 

C  WCU3D  IT  HIT  THE  OUTER  PLATE  IF  IT  WENT  PAR  ENOUGH? 

ELSE  IFfCHIP.GT.O.  .  AND .  UP-UPSPS  .  LS.  JOP  .  AND  .  CHIA.GE.  PI .  AND  .CHIA.  LT 
-RAD270)  THEN 

C  DCES  THIS  HAPPEN  INSIDE  THE  QUADRISPHERE? 

ITEST-5 

C  WOULD  IT  HIT  THE  INNER  PLATS  IF  IT  WENT  FAR  ENOUGH? 

ELSE  IF(CHIP.3E. 0. . AND . UI? . LE. UP+UPSPS .AND .CHIC. GE. PI  .AND.  CHIC.  LT 
*RAD2'0)  THEN 

C  DCES  THIS  HAPPEN  INSIDE  THE  QUADRISPHERE? 

ITEST-7 
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;test= : 

■X7.  =  -Z.\’’CZ 

■.'R  =  C3''SC-£A*Sa*CG 
RHC  =  Ar.ANG  ,  XH ,  :;r  1 
i.'rDi? 

:  '-VRITG:  6  ,  3CC1  .  N'CXX,!  .  RB/U  ,  PHI/PJVD  ,?.3/  'JI?  ,  RB/  '  UP-^UPEPS  )  ,RB/UP,RB/ 

C  -  'JP-JPEPS:  .RB/'JCP,R3/'JE.CHI?/RAB,GHIA/RAD,  CHIC/RAD,  ITSST 
3CC:  ECHMAT.  :  X  .  G -1  .  E3 . 3  .  E3 . 2 , 5?3 . 3  ,  3X  ,  3F3 . 2 , 1 A  ) 

RET’JRN 

ENE 

3G3RCUTINE  EWKM ( NE , 2ESM , EKST , EXENC , EBAR , SIG , ERAN , MAXX ) 

;.mpl:cit  real*3  (A-h^o-ei 
2IMENS2CN  2ESM(») 

SAVE 

maxx=:esm( 1 ) 

c 

C  BETERME.vrE  EBAR,  SIGMA 
NN=0 
EBAR=0 . 

SIG^O . 

ES’JM  =  0 
scsM=a . 

30  10  II-1,.VE 
^^^J=^^N-^IESM(  II ) 

EE  =  EXST*  ( II-l )  ■‘EKI.VC 
^S'JM=iESUM+EY*IESM(  II ) 

SQSM»SQSM+EY-*YY*ISSM(II) 

10  CCNTI.VUE 

ie;nn.ge.3)  then 
SBAR=ESUM/NN 

SIG* ( SCSM-EBAR«EBAR*NN) / (NN-l ) 
lElSIG.GT.O. )  SIG*SCRT(SIG) 

ENDI? 

MAXC-i- 1 
MAXX»ISSM(1) 

30  50  :«2,NE 

lEdESMd)  .GT.MAXX)  THEN 
MAXX-ISSMd) 

MAXCH-I 

ENDI? 

50  CONTINUE 

HMAX * FLOAT ( MAXX ) / 2 . 

CHl-0 . 

CH2-NE 
DEN-1 . 

30  60  I-1,MAXCH 

IF(HMAX.GT.?LOAT(IESM(I) ) . AND. HMAX.LE. FLOAT (lESM (1+1) ) )  THEN 
DEN-IESMd  +  l)-IESMd) 

IF(ABS(DEN) .LT.l.E-6)  DEN-1. 

CHl-I-l+(HMAX-FLOAT(IESMd)  )  )/DEN 
ENDI? 

60  CONTINUE 

DO  70  I-NE,MAX(MAXCa,2) ,-l 

IF(HMAX.GT.FLOATdESMd)  )  .  AND. HMAX .  LS  .FLOAT  ( IESM(  I-l )  )  )  THEN 
DEN-IESMd-l)-rESM(I) 

CH2-I-l-(HMAX-?LOAT(IESMd)  )  )/DEN 
ENDI? 

70  CONTINUE 

CCHl -EKST+CTl *SKINC 
CCH2-EXST+CT2*EXINC 
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iRAN-  '  d-iZ-Zl-il  ■  'SXINC 

?.ET'JR.V 

Z.'ID 

3:;3RC!j’Tr.VE  slips,  EXZ.RP,3ET,  the,  :3EL,3Q.^X?- hop,  RIA,R0A,-IR,3ET 
-?MAX , RMAX . RMIN . R.30 . 327  0 . ITSST , MTSJ ) 

IMPLICIT  PEAL’S  A-H,0-I1 

LCGI CAL  ;mX , CMN , QT? , CBT , QTC , CBC , CEX 

SAVT 

CATA  .'JCVLL  /O/ 

I  y  ■  .'fCALL  .  EQ  .  0  '  THEN 
?I=ACOS.-I. ) 

PI  32”! . 5*?I 
RTCD-I30 . /PI 
NCALL=I 
END  IE 

RLEL=2EL/RTCD 

RSET»BET/RTCD 

RTHE=THE/RTOD 

EXC=CQ/(2.*S?) 

CC=EXC/SXS/COS(RBET) **2 

CUKT-CC-TAN ( R3ET ) ’SIN ( RDEL ) + ( 1 . -CC ) *CCS ( REEL ) 

RR=R?/UURT 

?MAX*?I/2. 

lEICC.NE.l.  )  ?MAX--.ATAN^TAN(RBET)/(1.-CC)  ) 

UUPM-CC-TAN ( RBET ) "SIN ( PMAX ) + ( I . -CC ) *COS ( PMAX ) 

L’L'CM-CC^-TAN  ( RBET )  *SIN ( PMAX )  -  ( 1 .  -CC )  *COS  (  PMAX ) 

U 1 3 0 -CC- TAN ( RBET ) *SIN (  Pl-RTHE )  +  ( 1 . -CC) *COS ( PI-RTHE ) 

U2 7 0-CC-TAN ( RBET ) *SIN ( PI 32-RTHE )  +  ( 1 . -CC ) *COS ( PI  3 2-RTHE ) 

RMAX-RP/UUPM 

RMIN-RP/UUCM 

R130-RP/U180 

R270-RP/U270 

PMAX • PMAX* RTOD+THS 

IE(?MAX.LT.O.  )  PMAX-PMAX4-360. 

IE(?MAX.3T.36Q.)  PMAX-PMAX-360  . 

3ETN-  ( 1 . -CC ) «  SIN ( RDEL ) *TAN { RBET ) *COS ( RDEL ) 

3ETD-CC-t-  ( 1 .  -CC )  *COS  ( RDEL )  -TAN ( RBET )  *SIN  ( RDEL ) 

3ETA-0. 

IE ( 3ETD . NE . 0 . )  3ETA-ATAN( BETN/BETD ) *RTOD 
IE(BETA.GT.90. )  3ETA-BETA-360 . 

I?(3ETA.LT.-90. )  3ETA-3ETA+360 . 

ITEST  LOGIC 

I?(RMAX.GT.RMIN)  THEN 
PMlN-PMAX+180 . 

ELSE 

RT-RMAX 

RMAX-RMIN 

PT-PMAX 

PMAX-PMAX+180. 

RMIN-RT 

PMIN-PT 

ENDI? 

IE(?MAX.LT.O. )  PMAX-PMAX+360. 

IE(?MIN.LT.O. )  PMIN-PMIN+360. 

IE(PMAX.GT.360. )  PMAX-PMAX-360. 

IE(?MIN.GT. 360  . )  PMIN-PMIN-360  . 

QMX- PMAX .  GT  .  THE  .  AND  .  PMAX .  LS .  THE+DEL .  AND .  RMAX .  GE . ROP 
CMN- PMIN . GT . THE . AND . PMIN . LS . THE+DEL . AND . RMIN . LE . RIP 
I?(GMX.AND.CMN)  THEN 
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(  -i  ,  »  ,£.NC=  100  )  .'TRP  ,  (  IN(  X  )  ,;iAl  :  ,  I«1  ,NRP  ! 

C.\X.£  VEXP'NRP.XN.RAI 

r« 

RXP’-SAi  XAP'  3  .  X  ;  ^ 

£C?=RAiIA?'-i.X'.  ) 

RXAaRA. XAP; 3 ,1 ; ) 

3lCA»RAi  XAP  '  4  ,  X  )  ) 

•VKXrS  f  6  ,  '  ;  X  XX  .  ■  '  diameter  of  entrance  aperture  »  '  '  ,  ?7 . 3  ,  '  ■  CM  ’  '  ) 

-RCA-RXA 

•^XTE.'b,'  5X,''XNNER  AND  OUTER  RADIUS  OF  INSTRUMENT  »  '  ’  ,  2F7 . 3  , 

'  OM’ ' 1 ' :  RIP.ROP 

•^XTE( 5 , ' ' 3X . ' ■ INNER  AND  OUTER  RADIUS  AT  130  DEG  »  • ' , 2F7 . 3 

-  RAl  XAP' 3  ,  3  )  )  ,RA<  IAP(4, 3  n 

■^XTE.'o,  '  ;7X,  ’  'INNER  AND  CUTER  RADIUS  AT  EXIT  APP  »  '  '  ,  2P7 . 3 

-  '  ;  RA(XAP'3,-l)  )  ,RA(XAP(4,4)  ) 

C 

RXCN- I  X . /RIP- X . /RCP ) 

SXCN-ENCT*  ( RCP-t-RI?  )  *RXON 
VINN-SXCN/  ( 1  .  ■►RIP/RCP ) 

VCUT-VINN-SKCN 

RETURN 

100  CONTINUE 

STOP  ' NXTIN  * 

END 

SUBROUTINE  '/EX?  ( NRP  .  IN ,  RR ) 

IMPLICIT  REAL*3  (A-H,0-2) 

DIMENSION  IN(*),RR(*) 

DO  10  II-NRP,2.-1 
RR(IN(II) )-RR(II) 

DR=.{RR(II)-RR(II-l)  )/FLOAT(IN(II)-IN(Il-l)  ) 

DO  20  JJ»IN(II)-1,IN(II-1)+1,-1 
RR ( JJ ) “RR ( JJ+1 ) -DR 
20  CONTINUE 
10  CONTINUE 
RETURN 
END 
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on  no  oo  on 


;mX»QMX . AND . ?MAX . IT . ?MIN 
;mN- .NCT.QMX 

end:? 

:T?  =.;MX  .  dr  .  .^P  .  3T  .  RCP  .  OR .  RS .  32 .  ROP 
;3T=«  ;mN  .  OR  .  R?  .  IT  .  RIP  .  OR .  RS .  02 .  RIP 

;TC»THE.0T. :a0 . .AND.THE^DEL.OE. lao . .AND.R180 .32. ROA. and. R180.it. 
-RCP 

;3C-THZ.  IT.  lao  .  .AND.THE-t-DED.32. 180  .  .AND.R180  .OE.RIA.AND.R180 .3T. 

—  RT? 

02.'<  =  rHE.OT.  170.  .AND.THE+DEL.32.270.AND.R270.DT.ROP.AND.R270.CT.RI? 
IF'CTP.CR.QTC)  THEN 

0  HITS  TCP  PLATS  WITHIN  THE  *  DEL  INTERVAL 

:?(QTP.AND. (PMAX.L2.130. . CR . RiaO . 3T . ROA ) )  THEN 
ITSST=i2 

ELSE  :?(QTC)  THEN 
ITESTa^ 

ELS2 

:test»6 

ENDI? 

ELSE 

:tsst»i 

ENDI? 

I?(QBT.CR.QBC)  then 
C  HITS  3CTTCM  PLATE  WITHIN  INTERVAL 

I?;QBT.AND. (PMIN. L2. 180. .OR.R180.lt. RIA) )  THEN 
JTEST- 3 

ELSE  IF(QBC)  THEN 
JTEST-3 
ELSE 

JTSST»7 
ENDI? 

ELSE 

JTEST-1 
ENDI? 

I?(ITSST..NE.l.AND.3TEST.NE.l)  1TEST-MIN( ITEST  JTEST) 

I? ( ITEST . EQ . 1 )  ITEST- OTSST 
RETURN 
END 

SUBROUTINE  GEOMA ( VOUT . VINN .RIP, ROP , RIA , ROA , ENOT ) 

IMPLICIT  REAL*3  (A-H,0-Z) 

DIMENSION  IN(6S) ,XA(65),XT(6S) ,RA(65) , 

-t-''T(65)  .IAP(4,4)  ,?A(65) 

CHARACTER* 80  IHED 

RTOD-130./ACOS(-1. ) 

READ  INSTRUMENT  HEADER  CARD 
OPEN( 4 , PILE- ' TAPS4 . INS • ) 

READ(4, ' (A80) ' ,END-100)  IHED 
WRITS(6, ' (1H1//5X,A) ' )  IHED 

READ  LOGIC  CARD 

READ(4,2001,SND-100)  NX,.NT,NL,.NV,NIA,.NOA,NCA,IPRT 
101  ECRMAT ( 415, 2X, 311,15) 

READ  INPUT  APERTURE  CARD 
.READ  (4,*,  END- 100)  lAP 

READ  R  CARD  I.NPUT  AND  EXPAND 
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